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ABSTRACT 

We use high-quality Very Large Array (VLA) images of the Fanaroff & Riley Class 
I radio galaxy 3C31 at six frequencies in the range 1365 to 8440 MHz to explore the 
spatial scale and origin of the rotation measure (RM) fluctuations on the line of sight to 
O ' the radio source. We analyse the distribution of the degree of polarization to show that 

the large depolarization asymmetry between the North and South sides of the source 
C/3 , seen in earlier work largely disappears as the resolution is increased. We show that 

the depolarization seen at low resolution results primarily from unresolved gradients 
in a Faraday screen in front of the synchrotron-emitting plasma. We establish that 
the residual degree of polarization in the short-wavelength limit should follow a Burn 
law and we fit such a law to our data to estimate the residual depolarization at high 
resolution. We discuss how to interpret the structure function of RM fluctuations in 
the presence of a finite observing beam and how to address the effects of incomplete 
' sampling of RM distribution using a Monte Carlo approach. We infer that the observed 

RM variations over selected areas of 3C 31, and the small residual depolarization found 
Q^ ' at high resolution, are consistent with a power spectrum of magnetic fluctuations in 

' front of 3C31 whose power-law slope changes significantly on the scales sampled by 

our data. The power spectrum P{f) can only have the form expected for Kolmogorov 
' turbulence [P{f) oc /^^^^''] on scales <5kpc. On larger scales we find P{f) ^f "^'^- 

We briefly discuss the physical interpretation of these results. We also compare the 

■ global variations of RM across 3C 31 with the results of three-dimensional simulations 
rS , of the magnetic- field fluctuations in the surrounding magnetoionic medium. We infer 

■ that the RM variation across 3C31 is qualitatively as expected from relativistic-jet 

models of the brightness asymmetry wherein the apparently brighter jet is on the 
near side of the nucleus and is seen through less magnetoionic material than the 
fainter jet. We show that our data are inconsistent with observing 3C 31 through a 
spherically-symmetric magnetoionic medium, but that they are consistent with a field 
distribution that favors the plane perpendicular to the jet axis - probably because 
the radio source has evacuated a large cavity in the surrounding medium. We also 
apply our analysis techniques to the case of Hydra A, where the shape and the size 
of the cavities produced by the source in the surrounding medium are known from 
X-ray data. We emphasise that it is essential to account for the potential exclusion 
of magnetoionic material from a large volume containing the radio source when using 
the RM variations to derive statistical properties of the fluctuations in the foreground 
magnetic field. 

Key words: galaxies: jets - radio continuum:galaxies - magnetic fields - polarization 
- galaxies:ISM - X-rays: galaxies 
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The hot intracluster medium in rich clusters of galaxies 
is magnetised, with field strengths which are important 
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Figure 1. Grey scale at 5.5 arcsec FWHM resolution showing the 
main region s of 3C 31 in total intensity at 1.4 GHz as defined by 
iLaing et al] | |2008|. fig. lb). This image shows the radio source in 
equatorial co-ordinates, while the displays of depolarization and 
rotation measure below have been rotated anti-clockwise by 19.7 
deg to make the inner jet axis vertical. 



for the suppression of heat conduction, e ven if not dy- 
namically sig nificant (|Carilli fc Tavloj|2002l . and references 
therein). The fields are apparent via synchrotron radia- 
tion from cluster halo sources and Faraday rotation of po- 
larized emission from embedded and background sources 
(the subject of the present paper). At high spatial resolu- 
tion, the rotation of the E-vector position angle is accu- 
rately proportional to the square of the wavelength and 
depolarization with increasing wavelength is small, indi- 
cating that a foreground medium is responsible. Imaging 
of these foreground Faraday rotation measure (RM) vari- 
ations across radio sources in clusters has recently been 
used to estimate the power spectrum of RM and, by im- 
plication, of magnetic-field fiuctuatio ns, in turn leading to 
improved estimates of field strength (lEnfilin fc Vogt 20031: 
Voet fc Enf31inll200 3'.'2005V Murgia et al.ll2004l : lGovoni et al l 
20061 : 1 Guidetti et a l. 2008). RM fiuctuations are also seen 



against radio sources in sparser environmen t s such as galaxy 
group s (e.g. IPerlev. Bridle fc WiUid Il984l : iLaing fc Bridy 
1 19871 ). but the structure and strength of the magnetic fields 
in these systems have not been studied in detail. 

RM variations in clusters of galaxies have normally been 
modelled assumin g that the magnetic field is tangled either 
on a single scale (lBurnlll966 l: 'Felten"l996ll or over a r ange 
of scales (e.g. lEnfilin fc Vog t 2003; Murgi a et al]|2004l 'l and 
embedded in a smoothly varying, spherically symmetric den- 



sity distribution with parameters derived from X-ray obser- 
vations. This interpretation is complicated by three effects: 

(i) It is usually impossible to locate a radio galaxy along 
the line of sight within a cluster, and the path length over 
which Faraday rotation occurs is therefore not well deter- 
mined. It often has to be assumed that the sources are in 
the cluster mid-plane. 

(ii) The sources are inclined to the line of sight and may 
be extended on scales comparable with changes in external 
density, so the Faraday depth will vary systematically over 
the radio structure. The observed correlation between de- 
polarization and jet sidedness (in the sense that the side 
of the source containing the brighter jet depolarizes less 
rapidly with increasing wavelength) is most easily inter- 
preted in this way: if the nearer jet is brighter as a re- 
sult of Doppler boosting, then its emission is seen through 
less magnetoi onic materia l and therefore suffers less Fara- 
day rotation (jLain 3 1 19881 1. Although the general associa- 
tion between depolarization asymmetry and jet sidedness 
is well documented in both weak and powerful sources 
(|Morganti et al ] Il997l : iGarrington fc Conwavl Il99lh . there 
are as yet few sources that exhibit asymmetric depolar- 
ization for which it has been clearly established that the 
Faraday rotation is indeed due to a foreground medium (al- 
though there are no known counter-examples). The best- 
observed examples are all in central cluster galaxi e s wit h 
coolin g cores: CygnusA (jP reher. Carilli fc Perlcv' 19871}; 
M87 jOwen. Eilek fc Keel|[l990l: HvdraA (T aylor fc Pcrle^ 
Il993h and Hercules A l|Gizani fc L^I^Il999h . 

(iii) Radio sources must interact with their environments, 
leading to local deviations from spherical symmetry in 
the density distribution. There is direct evidence for the 
presence of evacuated cavities in the hot plasma coinci- 
dent with the radio lobes of sources in cl usters, as pre- 
dicted by models of r adio-s ource evolution (|Scheuej[l97i : 
iMcNamara fc NulsenI |2003). Magnetic fields in the sur- 
rounding medium may be important in stabilizing the cavi- 
ties against Kelvin-Helmholtz and Raleigh- Taylor instabih- 
ties (|Dursi fc Pfrommed [20081 ). 

In the present paper, we image and model the 
Faraday rotation distr ibution across the bright FRI 
iFanaroff fc RilevI Il97j ) class radio source 3C31. This 
source is identified with the nearby elliptical galaxy 
NGC383, w hich is the brightest member of a rich group 
of galaxies (|Arp| Il966l : IZwickv fc Kowall Il968l ). We there- 
fore probe a different environment from the Abell clusters 
studied by previous authors. Hot gas associated with the 
galaxy has been detected on both the group and galac- 
tic scales by X-ray im aging (jKomossa fc Bohringeil Il999l : 
iHardcastle et al.ll2002D . The nucleus of the source is located 
at the centre of the hot gas distribution, simplifying the ge- 
ometry, and th e orientation is well co nstrained, at least for 
the inner jets (|Laing fc Bridldl2002ah . 3C 31 shows a pro- 
nounced depolarization asymmetry in the expected sense 
l|Burchlll97 9l: IStrom et al. l ll983l : lAndernach et al.|[l993 '). Be- 
cause it is br ight, highly polarized, and well resolved in both 
dimensions (|Laing et al.l 2008. and references therein), the 
variation of the Faraday rotation measure across it can be 
studied in unusually rich detail. We can make a stringent 
test of the hypothesis of foreground rotation and evaluate 
the spatial statistics of RM over different parts of the source. 
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As the jet orientation is well constrained, we have a unique 
opportunity to compare the global changes in RM fluctu- 
ation amplitude across the structure with the predictions 
of three-dimensional models, including the effects of depar- 
tures from spherical symmetry. 3C 31 therefore provides an 
opportunity to examine the kiloparsec-scale structure of the 
magnetoionic medium in an environment that may be more 
characteristic of the FRI population as a whole than the 
clusters which have been the subject of many earlier stud- 
ies. 

We will show that our three-dimensional analysis of 
3C 31 is most consistent with the hypothesis that the ra- 
dio source has evacuated large cavities which are essentially 
devoid of thermal plasma, although existing X-ray obser- 
vations of the NGC 383 group are not deep enough to im- 
age such cavities if they are present. This motivated us to 
apply our techniques to the RM distribution of a second 
source. Hydra A, in which the cavities are well characterised 
I Wise et al.ll2007D . 

Section [2] of the present paper presents the main fea- 
tures of our data on the distributions of the Faraday rota- 
tion and depolarization of the polarized emission from 3C 31, 
and Section [3] gives an overview of our analysis and general 
assumptions. Sections |4] and [5] discuss the structure of the 
magnetic field in the foreground screen that is responsible 
for the Faraday rotation measure fluctuations. Section |4] de- 
scribes our analysis of the two-dimensional spatial statistics 
of the RM fluctuations. Section[5]explores three-dimensional 
models of the magnetoionic medium consistent with this 
two-dimensional analysis. Section [6] presents an analysis of 
the RM distribution in Hydra A. Section [7] summarises our 
conclusions. Some mathematical details are given in Appen- 
dices |A] and |B] and an analysis of the eflects of partially- 
resolved foreground rotation is described in Appendix [Cl 

Throughout this paper, we take N GC 383 to have a 
redshift of 0.0169 (the mean va lue fr omlSmith et al.ll2000|. 
Huchra. Vogelev fc GeUerl Il999l and [Oe Vaucouleurs et al.l 
199ll ). We assume a concordance cosmology with a Hub- 
ble Constant Ho = 70 km s"^ Mpc"\ = 0.7 and Q.m = 
0.3. At the redshift of NGC 383, this gives a linear scale 
of 0.344 kpc arcsec "^. For Hydra A, we take z = 0.0549 
jSmith et al.ll2004h and a linear scale of 1.067 kpc arcsec 



2 THE OBSERVED DEPOLARIZATION AND 
FARADAY ROTATION IN 3C 31 

T he observations and their reduction have been presented 
bv lLaing et al.l (|2008t ). along with discussions of spatial vari- 
ations in the radio spectrum and salient features of the total 
intensity and the apparent magnetic-field structure - all in- 
trinsic to the source. The principal large-scale features of 
the source discussed in that paper are shown schematically 
in Fig. [T] The data were obtained under NRAO observing 
proposals AF236 and AL405 at six frequencies - 1365, 1435, 
1485, 1636, 4985 and 8440 MHz - with angular resolutions of 
1.5 and 5.5 arcsec. The high quality and good sampling of 
the and angular domains in these data make them well 
suited for analysis of the statistics of the Faraday rotation 
and depolarization across 3C31. 

Magnetised thermal plasma along the line of sight 
causes Faraday rotation of the plane of polarization of lin- 



early polarized radiation, quantified by the rotation mea- 
sure, RM [equations (lAlfl - (|A3p ]. If the thermal plasma 
is mixed with the emitting material or if there are varia- 
tions of foreground Faraday rotation across the beam, the 
emission tends to depolarize with increasing wavelength. We 
will show that all of our observations are consistent with the 
expected variation of polarization with wavelength for an al- 
most completely resolved foreground screen. In Section [2. 11 
we establish some essential theoretical results which under- 
pin this analysis. We discuss the observed depolarization, 
the evidence that the Faraday rotating medium associated 
with 3C 31 is entirely in front of the radio source, and the 
variations of RM across the structure in Sections 12.21 12.31 
and 12.41 respectively. 

2.1 The effects of partial resolution for a 
foreground screen 

The observing beam modifies the derived spatial statistics 
of RM. In practice, we derive E-vector position angle from 
Q and U images convolved with the beam and use the po- 
sition angles at multiple frequencies to calculate RM. The 
beam has the effect of suppressing high spatial frequencies 
but for wavelengths close to (and consequently low depo- 
larization), this effect can be evaluated simply. The argu- 
ment, whi ch is fundamen tal to the analysis in Section [l] is 
as follows (|Tribblelll991al ). For a complex polarization 

p(A) = |p(A)|exp[2ix(A")] 

defined to have zero position angle for A = 0, the RM at any 
wavelength is given by 



_ 1 /ainpN 

In the case of pure foreground rotation, |p| is independent 
of A. If, in addition, it does not vary over a beam W{r±) 
of unit area, i.e. p(rx,A) = po, where po is a constant and 
rx is a vector in the plane of the sky, then the measured 
polarization is 

p(rx,A) = |po|y"w^(rx-r;)exp[2iRM(r'x)A']d2r'x 

The observed RM is derived by differentiating under the 
integral sign and taking the zero- wavelength limit: 

RM(rx) = y W(rx - r'x)RM(r'x)d^r'x (1) 

i.e. the measured RM distribution is closely approximated 
by the convolution of the true RM distribution with the ob- 
serving beam. We use this result to include the observing 
beam explicitly in our analysis, verifying post hoc that the 
approximation is valid. 

This approach also makes it easy to see why A^ ro- 
tation persists even in the presence of small but signif- 
icant depolarization. Consider a circular Gaussian beam 
M^(rx) = exp(-ri/2a^)/27rcr^ To first order 

RM(rx) ^ RM(0) + rx-VRM 
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Figure 2. Depolarization and rotation measure at a resolution of 5.5 arcsec FWHM. (a) Burn law k in rad^ m"**, derived from a weighted 
least-squares fit to the relation ln[p(A)] = ln[p(0)] — kX^ between 4985 and 1365 MHz. The colour range is from —400 to 1000 rad'^ m~^. 
(b) RM from a fit to the four frequencies between 1636 and 1365 MHz, primarily to show the RM in the North tail and spur. The range 
is —200 to — 33radm~'^. (c) as (b), but using 5 frequencies between 4985 and 1365 MHz. The errors are smaller than in panel (b), but 
fewer points in the North of the source have sufficient signal to determine the RM. The rectangles marked 'NP2' and 'SP3' denote the 
regions for the RM histograms in Fig.lSfa) and (e) and the structure-function analysis of Section 14.41 The crosses mark the positions for 
the plots in Figs |4] and |6] 



where VRM is also evaluated at rx = 0. The integral for 
p(A) can then be evaluated analytically to give: 

^ |p|exp[2i(RM(0) + rx-VRM))A^] 

X exp[-2|VRM|V^A''] (2) 

The first exponential term corresponds to rotation and 
the second to depolarization. Thus if the RM varies by a 
small amount across the beam, then we expect A^ rotation 
even though there is significant (albeit small) depolariza- 
tion. The wavelength depe ndenc e of de polarization has the 
functional form derived bv lBurnI (|l966l ): 

p(A) =p(0)exp(-fcA*) (3) 

with k = 2| VRM 1^0"^, but no te that this wil l not in general 
persist at longer wavelengths l|Tribblelll991al . see Section [431 
and Appendix [C|) . 



2.2 Depolarization 

The effect of depolarization between two wavelengths Ai 
and A2 is usually expressed in terms of the ratio DP^J = 
p{\2)/p{\i), where p(A) is the degree of polarization. Our 
six-frequency data set is capable of providing more sophis- 
ticated information about the depolarization, but to make 
effective use of it we need a suitable fitting func tion. As we 
have shown in Section 12.11 a Burn law variation |Burnlll966l : 
equations and Q], is expected for partially resolved fore- 
ground Faraday rotation in the short-wavelength limit. The 
data at the vast majority of individual pixels are fit to within 
the errors by this relation and we show images of k at 5.5 
and 1.5-arcsec resolution in Figs. [2ja) and|3ja)Q Errors in 

^ The wavelength dependences at a small minority of points in 
the 5.5-arcsec images of the Southern part of the source show 




500 1000 -150 -100 -50 -150 -100 -50 



S 



a) Burn aw k 





1 \ \ r 



1 \ \ \ \ \ r 



fc) RM 



NP1| 



SP1 



J \ \ L 



J L 



80 60 40 20 -20 -40 
ARC SEC 



80 60 40 20 -20 -40 
ARC SEC 



80 60 40 20 -20 -40 
ARC SEC 



Figure 3. Depolarization and rotation measure at a resolution of 1.5 arcsec FWHM. (a) Burn law k in rad^ m derived from a weighted 
least-squares fit to the relation ln[p(A)] = ln[p{0)] — fcA'* for six frequencies between 8440 and 1365 MHz. The colour range is from —400 
to 1000 rad^ m as in Fig. (Sfa). (b) RM determined by a fit to position-angle images at 6 frequencies between 1365 and 8440 MHz 
using a version of the AlPS task RM modified by G.B. Taylor. The colour range is —180 to +20radm~^. The boxes NPl, SPl and SP2 
show the areas used fo r the RM histograms in Fig. [8] and for the structure-function analysis in Section 14.41 (c) As (b), but using the 
PACERMAN algorithm of lDolag. Vogt fc EntilinI ||2005| ). The crosses mark the positions for the plots in Figs [5] and [7] 



k due to rotation across the observing bands arejwgligible 
for th e bandwidths used for these observations (jLaing et al.l 
|2008| . table 1) and the observed Faraday rotations (see Sec- 
tion [23J. 

The well-established depolarization asymm etry between 
the N orth and South sides of 3C 31 (Burch 1979: IStrom et all 
is still present at a resolution of 5.5 arcsec, although 
reduced in amplitude from that in the lower resolution 



significant deviations from a Burn law. We return to this and 
related issues in Section 14.51 and Appendix [Cl 



data. Fig. [IJa) shows k derived from a fit to 5-frequency 
data between 1365 and 4985 MHz at this resolution. Exam- 
ple plots of p against A** are shown in Fig. |4] in order to 
demonstrate that the Burn law gives adequate fits for this 
frequency range and resolution. Significant depolarization 
(fc > 0) is found at most individual points, but a minor- 
ity have fc < due to noise or the effects of Faraday rota- 
tion on a non-uniform distribution of intrinsic polarization. 
The mean values for the North and South of the source 
are (fc) — SSrad^m"* and (fc) = 364rad^m~'', respectively, 
corresponding to DPl^JjJ" ^ 0.8 and 0.4, compared with ~1 
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Figure 4. Plots of degree of polarization, p (log scale), against 
A* at the positions indicated by crosses in Fig. |2] and quoted in 
the individual panels. The resolution is 5.5arcsec FWHM. Burn 
law fits [equation ([3j] are also plotted. 



Figure 6. Plots of E-vector position angle x against at the 
positions indicated in Fig. |2] The resolution is 5.5arcsec FWHM. 
Fits to x(A2) = x(0) + RMA2 are shown, and the values of RM 
are given in the individual panels. The ranges in angle are 300° 
in panels (a) - (b) and 600° elsewhere. 
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Figure 5. Plots of degree of polarization, p (log scale), against 
A* at the positions indicated by crosses in Fig. |3] and quoted in 
the individual panels. The resolution is l.Sarcsec. Burn law fits 
[equation ((Sjl] are also plotted. Of the fits in this figure, only (a) 
and (c) show significant depolarization. 
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Figure 7. Plots of E-vector position angle x against A^ at the 
positions indicated by crosses in Fig. |3] Fits to x(-^^) = X(0) + 
RMA^ are shown, and the values of RM are given in the individual 
panels. The resolution is 1.5 arcsec FWHM and the position-angle 
range is 300 deg. 
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and 0.1 - 0.3 at re solution of 24arcsec x 47arcsec FWHM 
jStrom et al.lll983l ). 

Comparison between Figs [2ja) and |3ja) shows further 
that there is much less depolarization at 1.5-arcsec resolu- 
tion. The data at this resolution are noisy and the fitted 
values of k are consistent with zero at the majority of indi- 
vidual points (many points therefore have A: < 0). Of the sky 
positions shown in Fig. O only (a) and (e) show significant 
depolarization, for example. Positions with fc > tend to 
occur in coherent clumps (the two most prominent are la- 
belled DPI and DP2 in Fig. [3^), indicating that the residual 
depolarization in the South is still significant, however. Sim- 
ulation of our image analysis and fitting procedures shows 
that the distribution of errors in k is position-dependent and 
significantly asymmetrical over much of the source at 1.5- 
arcsec resolution, leading to biases in spatial averages. For 
quantitative comparisons, we therefore use only data with 
p > Arjp at all frequencies, in which case we estimate sys- 
tematic offsets in the mean values to be negligible compared 
with the quoted random errors. Averages over these posi- 
tions confirm not only the decrease in depolarization with 
increasing resolution, but also the persistence of a residual 
difference between the two jets. In a high signal-to-noise re- 
gion within 60arcsec of the nucleus in the North jet, we 
find {k) = 3 ± 7 rad^ m~* compared with 32 rad^ m~^ at 
5.5 arcsec resolution. Similarly, for the inner 60 arcsec of the 
South jet we find (k) = 55 ± 7rad^m~* (1.5 arcsec) and 
166rad^m~* (5.5 arcsec). 

Apart from the large-scale North-South asymmetry, 
there is no compelling evidence for any detailed correlation 
between depolarization and source structure (as might be 
expected if the effect occurred in a thin surface layer just in 
front of the emitting material). 

2.3 Faraday rotation: images and evidence for 
foreground rotation 

RM images at 5.5-arcsec resolution made using a version 
of the AlPS task RM modified by G.B. Taylor are shown 
in Fig. \^h) for the 1365-1636 MHz data alone and in 
Fig-El^c) for all five frequencies between 1365 and 4985 MHz. 
The former is noisier in the region of overlap, but covers a 
wider area: the two images are consistent with each other 
where they overlap. Example fits of E- vector position angle 
X against are shown in Fig. [S] 

There is a clear asymmetry in RM fiuctuation ampli- 
tude, which is much greater in the South. In the outer North 
jet, spur and tail, the low depolarization (Fig. [5^) and ro- 
tation together indicate mostly resolved foreground Faraday 
rotation. For the parts of the North tail and spur which show 
significant polarized signal at 1365 to 1636 MHz but which 
are not detected at 4985 MHz, the depolarization found at 
lowe r resolution betw een these frequencies is also negligi- 
ble (jStrom et al.lll98^ ). Large-scale coherent RM variations 
(«10 radm~^ on scales from «10 arcsec up to «larcmin) 
are observed in this part of the source [Figs [2][b) and (c)]. 
These are large enough to account for the depolarization of 
0.3-0.4 observed between 1.4 and 0.6 GHz with a beam of 
29 X 55 arcsec (|Strom et al.lll983l '). In the South spur and 
tail, despite the large RM gradients and significant depolar- 
ization, the rotation remains accurately proportional to A^ 
except at a few positions. These are usually associated with 




RM / rad m-^ 

Figure 8. Histograms of the RM distributions within the regions 
marked NP2, SP3 (Fig.©, NPl, SPl and SP2 (Fig.©, (a) NP2. 
The RM distributions from 5- and 4- frequency fits at 5.5 arcsec 
FWHM are shown by full and dotted lines, respectively, (b) — (d) 
NPl, SPl and SP2. The distributions at 1.5 arcsec FWHM from 
the AlPS and pacerman fitting algorithms are shown by the full 
and dotted lines, respectively, (e) SP3, showing the RM distribu- 
tion from a 5-frequency fit at 5.5 arcsec FWHM. The histograms 
have been scaled independently to approximately the same max- 
imum value and the horizontal error bars show the mean and 
maximum rms errors derived from the RM fits. 



abrupt changes in RM on scales comparable with the observ- 
ing beam and high depolarization. We show in Appendix [Cl 
that these effects are quantitatively consistent with Faraday 
rotation by a partially resolved foreground screen. 

Images of rotation measure in the jets at 1.5-arcsec 
resolution using all 6 frequencies are shown in Fig. O We 
used two different fitting programs: the modi fied RM task 
[Fig- Elb)] and the pacerman algorithm of iDolag etliD 
(i2009 . Fig. [3J;) . In the former case, we required that all E- 
vector position angles at a given point had rms error < 10° , 
but still found that a few points had wildly discrepant RM 
values due to erroneous resolution of mi ambiguities. Un- 
like the effect noted at 5.5-arcsec resolution (where the RM 
anomalies are extended over several beamwidths), the er- 
rant values occurred at isolated pixels and cannot represent 
real RM structure. They were therefore edited out. The pac- 
erman algorithm solves for the nn ambiguities in regions of 
high signal-to-noise ratio and uses this information to resolve 
them in adjoining, fainter areas. Where RM solutions were 
found by both algorithms, they were in excellent agreement, 
but pacerman gave credible results for many more points 
with low signal-to-noise ratio, especially in the South. 

At 1.5-arcsec resolution, there are no significant devi- 
ations from the relation oc A^ over a range of rotation 
Ax up to ~300° (fits at representative points are shown 
in Fig. (Tjl. Taken with the observation that the depolar- 
ization is very slight over our frequency range even in the 
South (Section 12. 2|) . this implies that most, if not all of the 
Faraday rotation of the jet emission must be produced by 
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Figure 9. Profiles along an axis defined by the inner jets in 3C31 (PA 19.7°). Positive distances are in the direction of the North jet 
and the vertical dotted lines show the position of the nucleus. Quantities are evaluated for boxes of length 65 arcsec along the axis for 
the 5.5 arcsec FWHM data (panels a, c and e) and 20 arcsec for the 1.5 arcsec FWHM data (panels b, d, f), using the pacerman fitting 
algorithm for the RM in the latter case. The boxes are wide enough to include all unblanked points, (a) and (b) Mean RM, (RM). The 
horizontal dotted line shows our adopted value of the foreground rotation measure, RMq. (c) and (d) ctj^Mi determined with respect to the 
local mean in the box. The values of ctrm have been corrected to first order for the offset due to noise on the fit. (e) and (f) Mean value 
of the Burn law parameter k, with statistical uncertainties estimated from the least-squares fits. Positive values imply depolarization 
with increasing wavelength and k = 100 rad^ m~'* corresponds to a depolarization ratio of 0.79 at 1365 MHz. 



foreground material and that most of the depolarization ob- 
served at lower resolutions results from RM fluctuations on 
scales larger than 1.5 arcsec. The residual depolarization at 
1.5 arcsec resolution is probably caused by RM fluctuations 
on even smaller scales. We cannot establish this conclusively 
without higher-resolution data, but we show consistency of 



the observed wavelength dependence and magnitude of the 
depolarization with models of the spatial variation of fore- 
ground RM in Sections 12.11 and 14.51 respectively. The most 
extreme depolarization (e.g. regions DPI and DP2) tend to 
be associated with large RM gradients, as expected for fore- 
ground models [Figs[3ja) and (c)]. 
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2.4 Variation of Faraday rotation and 
depolarization along the jets 

The RM distributions for regions NPl, NP2, SPl, SP2 and 
SP3 are plotted as histograms in Fig. [S] The results for 4- 
and 5-frequency fits are in good agreement for region NP2 
[Fig. IHa)], although the former are inevitably noisier. The 
MPS and PACERMAN algorithms give very similar results for 
NPl and SPl but they disagree significantly in SP2, where 
the former fits many fewer points [Fig. IHd)]. We use the 
1.5-arcsec pacerman images and the 5.5-arcsec, 5-frequency 
fits for quantitative analysis. Note that the fitting errors 
(indicated by the error bars in Fig. [8]) are generally much 
smaller than the widths of the histograms. 

In the North, the RM distributions are sharply peaked 
with a total range from —100 to — 40radm~'^. The inner N 
jet shows a similar range. In the South, the RM range is 
much larger (—200 to +30radm~^), and the distributions 
appear more uniform. 

Fig. |9] shows profiles of the mean rotation measure, 
{RM), the rms with respect to this local mean, ctrm and the 
parameter k in the Burn depolarization law along the inner 
jet axis (position angle 19.7 deg). The 5.5- and 1.5-arcsec 
data were averaged over boxes wide enough to contain all 
unblanked pixels and of lengths 65 and 20 arcsec along the 
axis, respectively. In order to avoid systematic errors, only 
points with p > 4crp were used to generate the profiles of k. 
We made a first-order correction to the rms RM, o-RMraw, 
by subtracting the fitting error a^t in quadrature to give 
o"RM = (fRMraw " ''"It)^''^! but the cffect of this correction is 
in any case very small. 

There is a systematic gradient in {RM} along the North 
jet, but the values vary erratically in the South, ctrm in- 
creases from North to South from a minimum of ~5 rad m~^ 
at 5.5 arcsec resolution (300 arcsec North of the nucleus) to 
a maximum of 30-40 rad m"'^ (depending on resolution) at 
200 arcsec to the South. There is only evidence for a decrease 
in the last box 400 arcsec from the nucleus in the South, 
where our data are limited in signal-to-noise and restricted 
to four observing frequencies. The profile of (k) at 5.5-arcsec 
resolution [Fig.[9je)] shows that depolarization also increases 
monotonically from North to South across the nucleus, but 
there is evidence for a decrease in the South at larger dis- 
tances. At 1.5 arcsec resolution, depolarization is significant 
only in the South jet. The differences between the two res- 
olutions are as expected for partially resolved RM fluctua- 
tions: at higher resolution there is less depolarization and 
lower rms but variations between box means are larger. The 
increase in RM fluctuation amplitude going from North to 
South across the nucleus is striking. 

The very small variation in {RM) in the North of the 
source suggests a Galactic origin. Estimates of the Galac- 
tic foreground RM derived from observ ations of nearby 
sourc es range from —64 to — 39radm~^ (|Andernach et al.l 
Il992|) and the val u e pre dicted by the smoothed models 
of iDineen fc ColesI (120051 ) is -49.3 rad m"^ both consis- 
tent with this idea. The best estimate of the foreground 
RM from our own data comes from the region more than 
200 arcsec North of the nucleus of 3C 31, where fluctuations 
are smallest [Fig. HJa)]. In what follows, we adopt the un- 
weighted mean RM from the 5-frequency fits over this re- 
gion, RMo = — 56.9radm~^, as the Galactic contribution. 



Note, however, that we cannot exclude a significant large- 
scale component of RM local to this part of the source (Sec- 
tion I4.4|l . 

2.5 Limits on internal Faraday rotation 

Our limits on depolarization and deviations from rota- 
tion, within 60 arcsec of the nucleus in the North jet can 
potentially constrain the internal density of thermal elec- 
trons. For a line of sight intersecting the jet axis 28 arcsec 
from the nucleus in projectio n (the maximum exte nt of 
the conservation- law model of iLaing fc Bridie l2002bl ). our 
kinematic model gives a path length through the jet of 
7.3 kpc an d an on-axis equipartition magnetic field of 1.5 nT 
(|Laing fc Br idle 2002a.) . We adopt these as fiducial numbers, 
noting that there will be variations in path length and field 
along and transverse to the jet axis and that the component 
of the magnetic field along the line of sight will be reduced 

due to projection. 

We use the formulae given by iBurnI l| 19661 ). taking a 
conservative upper limit of {k) < 20rad^ m"'* (equivalent to 
a depolarization of >0.95 at our lowest observing frequency 
of 1365 MHz). For a uniform magnetic field B along the line 
of sight and a path length L, this corresponds to a limit on 
the internal density of 

n/m"^ < 120(B/1.5nT)~^(L/7.3kpc)"^ (4) 

The lack of deviation from a law requires that no more 
than ~45° of the total rotation is due to thermal matter 
distributed uniformly within the jet, giving n < 370 m~^ 
for the fiducial parameters and a uniform field. If, on the 
other hand, we assume that the field component along the 
line of sight has a Gaussian distribution with rms B and A'^ 
reversals through the jet, we find 

n/m"^ <40(B/1.5nT)"^(L/7.3kpc)"^A'^''^ (5) 

Our modelling of the magnetic-field structure in the jets 
iLaing fc Bridle! l2002ah implies that the magnetic field is 
dominated by the toroidal component, with a smaller lon- 
gitudinal contribution. There are two extreme possibilities 
iLaing. Bridle fc Canvinll2006D : 

(i) The dominant toroidal component is vector-ordered, 
but the longitudinal component has many reversals. We 
would then expect systematic variations in Faraday rotation 
transverse to the jet axis, which we do not observe. 

(ii) The field is entirely disordered, but anisotropic. 

In either model, there would be systematic (and so far unob- 
served) variations of depolarization across the jets, together 
with a gradient in RM if the toroidal component is ordered. 
A pure toroidal field has no component along the line of 
sight on the projected jet axis, so we would expect to see 
little Faraday rotation and depolarization there. The density 
limits for there more realistic field configurations would be 
a few times larger than those given by equations Q and ((5)1 
for cases (i) and (ii), respectively. The number of reversals, 
A'', in the fully disordered case could be very large. 

Even in the case of a fully ordered field, these lim- 
its are fully consistent with the densities inferred from our 
conservation-law analysis on the assumption that the ther- 
mal electrons originate entirely from entrained plasma with 
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normal cosmic abundances: n ~ 0.7 m ^, also at a pro- 
jected distance of 28 arcsec from the nucleus (|Laing fc Bridle! 
l2002bl ). 



3 OVERVIEW OF THE ANALYSIS 
3.1 Description 

We seek to interpret the depolarization and Faraday ro- 
tation results in terms of the magnetic field structure in 
the Faraday-rotating medium, first by deriving the two- 
dimensional statistics of the magnetic-field fiuctuations that 
produce the observed Faraday rotation (Section|4)) and then 
by constructing models of three-dimensional distributions 
of gas and magnetic field which are consistent with these 
statistics (Section [Sjl. 

Our analysis proceeds in three stages: 

(i) The two-dimensional fluctuations of Faraday rotation 
measure on scales larger than the observing beamwidth can 
be quantified in several ways. We will justify use of the RM 
structure function (Section 14.21) as the principal statistic. 
This restricts us to regions small enough that lateral vari- 
ations of the foreground density and rms field strength can 
be neglected. Our approach is to guess a form for the mag- 
netic power spectrum, derive the corresponding RM struc- 
ture function (modified by the observing beam) and compare 
it with the data, estimating errors for the poorly-sampled 
large separations by making a set of artificial RM images 
with the theoretical power spectrum sampled on the same 
grid as our data. 

(ii) Residual depolarization can be used to estimate the 
variations of F araday rotation across the beam (Burn 19 66: 
lTribble|[l991ah . thus extending the range of spatial scales we 
can study. We evaluate the depolarization expected from dif- 
ferent power spectra numerically, by making finely-sampled 
realisations of RM images, which we then use to generate 
artificial polarization maps at the observing resolution. 

(iii) For a large-scale foreground medium, the varia- 
tions of RM fluctuation amplitude and depolarization 
with position depend on t h e orientation of the source 
l|Garrington fc Conwavl [l99ll : iTribbld [l99? ) . Knowing the 
density distribution from the X-ray observations and the 
geometry, the strength of the magnetic fleld can be esti- 
mated by fltting to the observed distributions using an- 
alytical single-scale models (e.g. [Felton 1996), appropri- 
ately weighted estim ates of the au t ocorr e lation function 
or p ower spectrum (|EnfiUn fc Vod l2003l: IVogt fc Eni3Unl 
20051) or three-dimensional simulations~(|Murgia et al.ll2004l : 
Govoni et all l2006l : iGuidetti et"ai1 12008| ). We use simula- 



tions, incorporating constraints on the magnetic-field power 
spectrum from (i) and (ii). Unique features of our analysis 
are that the orientation of 3C 31 (at least close to the nu- 
cleus) is well determined and we have enough data to test 
whether the Faraday-rotating medium is spherically sym- 
metric with a radial profile related to the density of the 
surrounding hot plasma. 

The advantage of the approach using two-dimensional RM 
and depolarization analyses [steps (i) and (ii)] is that when 
making artificial RM images, we can use large computational 
grids to sample a full range of spatial scales accurately and 



efficiently, reserving three-dimensional simulations for stud- 
ies which include global variations of density and magnetic- 
fleld strength. We can also use the Hankel transform relation 
fSection l4.2.1|l to predict the structure function expected for 
a given power spectrum, allowing explicitly for the effects of 
the observing beam. 



3.2 Simplifying assumptions 

We make a number of simplifying assumptions, as follows. 

(i) Faraday rotation is produced entirely by a foreground 
screen. We have argued in Section 12.31 that this is a very 
good approximation for 3C 31. 

(ii) The magnetic fleld is a Gaussian random variable, 
and its spatial distribution can therefore be described by 
the power spectrum or its Fourier transform, the autocor- 
relation function. This may well be incorrect if the field 
is intermittent, e.g. in the form of flux tubes or filaments 
(|Eilek fc Owen|[2002l . and references therein), but our data 
are too sparse to allow us to determine higher-order corre- 
lations. 

(iii) The fleld is isotropic, in the sense that it has no 
preferred direction when averaged over a sufflciently large 
volume. Our data show no evidence for a preferred direc- 
tion in the RM distribution on scales < 100 arcsec, but we 
note that this may not al ways be the case: for example the 
RM distributions in M 84 (|Laing fc BridldllQaTl ) and 3C 465 
(|Eilek fc Owen|[200^ ') appear to have banded structures. We 
do see RM gradients on very large scales in 3C 31, but these 
are consistent with sampling of an isotropic magnetic power 
spectrum with significant amplitude on these scales together 
with spatial variations in density. 

(iv) The amplitude of the magnetic-field power spectrum 
varies with the thermal gas density in the foreground screen, 
which is in turn a smooth function of position. The form of 
the power spectrum, in contrast, is independent of position. 
We test this assumption by evaluating the spatial statistics 
in different regions of the source. 

The importance of these assumptions is that the power spec- 
tra of the magnetic field and RM are then proportional 
lEnfilin fc Vogtll2003l 'l . We work in terms of spatial frequency 
f = (^i/i/i/z) with fx = {fx,fy) in the plane of the sky. 
In the isotropic case, we can define the three-dimensional 
magnetic power spectrum w{f) such that w{f)dfxd[ydfz is 
the power in a volume df^dfydf^ of frequency space|j 

Similarly, the two-dimensional RM power spectrum 
C{f±) is defined such that C{f±)dfxdfy is the RM power 
for an element dfxdfy. For a foreground screen of uniform 
depth and constant electron density, w{f) and C{f±) are 
simply related [equation ()A6p ]. If the electron density and 
the normalization of w{f) vary significantly along the line 
of sight but not perpendicular to it, the functional form s 
of w{f) and C{f±) are still identical (|Enfilin fc Vogtll2003l ). 
Thus we may use the RM structure function evaluated over 
limited areas to determine the form of w{f ), but estimating 



^ Note that some authors use the one-dimensional form 
47r/^u)(/). For the three-dimensional form, as used in this pa- 
per, a Kolmogorov power spectrum is w{f) oc f~^^^''^ whereas 
the one-dimensional form is iiTf^u){f) oc f~^^^ ■ 
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its normalization as a function of position requires a three- 
dimensional analysis. 



4 TWO-DIMENSIONAL ANALYSIS 
4.1 Choice of analysis technique 

Our goal is to characterise the RM over areas of the radio 
source where the parameters which define its spatial statis- 
tics do not vary significantly. Descriptions in terms of spa- 
tial or frequency variables are equivalent and several dif- 
ferent techniques have been used in the lite rature. One ap- 
proach is to work directly in frequency space (lEnfilin fc Vogtl 
I2OO3I : IVogt fc EnBhnI |2003| ) to determine the power spec- 
trum. Real-space alternatives are: 

(i) the autocorrelation function 
C(rx) = {RM(ri + rY)RM(rl)> 

lEnfilin fc Vogtll2003l : IVogt fc Enfiiinll2003l ) and 

(ii) the related structure function defined by: 

S{r^) = ([RM(ri + r^) - RM(rY)]'> 



(e.g. ISimonetti. Cordes fc Spangleil Il984l : iLeahvl 1 19871 : 
iMinter fc SpangleijligQef T 



Here, the vectors and are again in the plane of 
the sky and {) denotes an average over r^, defined by the 
area of the radio source. For a stationary random process, 
the structure and autocorrelation functions are related by 
S{r±) = 2[C(0) — C{r±)], although in practice, this relation 
will be approached only if the largest scale of fluctuations is 
significantly smaller than the size of the measurement area. 

We work in real space and use the structure function 
rather than the autocorrelation function for the following 
reasons: 

(i) Most of the useful properties of the autocorrelation 
function only appear when C{r±) drops to at some outer 
scale small compared with the sampled area. If this is not the 
case (as in 3C31), the form of the function can vary wildly, 
typically oscillating with large amplitude at large separa- 
tions. 

(ii) A related point is that sensible use of the autocorre- 
lation function requires an accurate knowledge of the zero 
level. This is uncertain for 3C31: we have argued that the 
mean RM in the North tail is primarily Galactic, but there 
could still be significant contributions from local material. 
The zero-level for a particular region cannot be determined 
merely by averaging in the presence of large-scale fiuctua- 
tions. 

(iii) This problem is made worse by the irregular shapes 
of the sampling regions in 3C31. This leads to particular 
problems with the frequency -space approach, as the corre- 
sponding window functions (jEnfilin fc Vogtl [20031 ) are very 
complicated. 

(iv) By comparison, the structure function S{r±) is inde- 
pendent of the mean level and (to first order) of structure 
on scales larger than the measurement area. 

S{r±) is therefore relatively robust and can use all of the 
available evidence without giving systematic biases at the 
largest separations. 
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Figure 10. Full lines: plots of the structure functions correspond- 
ing to power spectra C(/x) oc /j^' (for p = 2.6, 3.0 and 3.4) over 
an infinite frequency range, from equation ||B1| |. Dotted: structure 
functions for the same power spectra, but for images convolved 
with a Gaussian of FWHM 5 pixels [equation ||B2| |] . The crosses 
are the results of numerical integration. The agreement between 
analytical and numerical results is very good except for the case 
of a flat power spectrum without convolution, where there are 
small but significant errors at small separations (this case is not 
relevant to the observations, however). 



4.2 Theoretical structure functions 

4-2.1 Hankel transform relations 

For isotropic fluctuations, the RM power spectrum C{f±) is 
the Hankel tr ansform of C(r±) . Using the Fourier transform 
convention of iBracewelll l|2000l ). 



C{ri 



2n I C(/i)/xJo(2^/xri)d/_L 

10 



(6) 



where C{f±) depends only on the scalar spatial frequency 
/x . In the short- wavelength limit for a beam of FWHM — 
2a{ln2/Tv)^^^ , the relation becomes 



C(rx) = 27r / C'(,fx)/±Jo(27r/xr_L) 



X exp( — 27r(j fi_)df± 



(7) 
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4-2.2 Power-law power spectra and the effects of the 
observing beam, 

In what follows, we use finite frequency intervals and there- 
fore need to evaluate the Hankel transform integrals [equa- 
tions (|6]) and ((7|)] numerically. As a test of our procedure, 
Fig. [To] shows a comparison of analytical and numerical es- 
timates for structure functions corresponding to power-law 
C{f±) for an infinite frequency range, both before and after 
Gaussian convolution of the image. The analytical expres- 
sions are given in equations (|B1[I and (|B2[l . 

The effects of convolution on the observed structure 
function are very noticeable, especially for flat power spec- 
tra. For example, the structure function for q — 2.6 deviates 
significantly from its asymptotic form even for separations 
r± « 5xFWHM (Fig. [10)). After convolution, the structure 



function approaches S{r±) oc rj_ for small r± (Section [B2] 



4.3 Realisations of the RM distribution 

4.3.1 Method 

Errors for the structure function of poorly sampled data are 
impossible to evaluate analytically: they are correlated be- 
tween adjacent position bins and depend on the details of 
the window function, particularly for large separations. For 
these reasons, we choose to generate multiple realisations of 
Gaussian random RM fields with given power spectra on the 
observed grids. The dispersion in the structure functions cal- 
culated for these realisations gives an estimate of the errors 
and allows us to determine the range of scales over which 
the observed structure functions are reliable. Our method 
for generating a realisation of a Gaussian random variable 
with a particular power spectrum is as follows. 

(i) Make an array in the frequency plane whose elements 
are complex numbers whose real and imaginary parts are 
random numbers selected from a Gaussian distribution with 
zero mean and unit variance. 

(ii) Multiply the elements of this array by the square root 
of the power spectrum, [C{f±)]^^'^. 

(iii) Optionally, multiply by the Fourier transform of the 
convolving beam. 

(iv) Perform a two-dimensional inverse Fourier transform 
to the image plane. 

(v) Evaluate the structure function. 

This was implemented using a EFT routine on a grid of size 
n = 2™, setting implicit limits 1/n < f± < 1/2 to the max- 
imum and minimum frequencies which can be sampled in 
either dimension (the upper limit is the Nyquist frequency) . 
In practice, we used n — 8192 to ensure that the difference 
between the structure functions for minimum frequencies of 
and 1/n is negligible over the range of spatial separations 
sampled by our data. 
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Figure 11. Comparison of structure functions derived by simu- 
lation and Hankel transform for an observing beam of 1 . 5 arcsec 
FWHM. (a) Broken power-law with indices 5high = 11/3, giow 
= 2.32 and a break frequency /i, = 0.062 arcsec"^ [equation ISllI . 
(b) Power law with index q = 2.39 and a high-frequency cut-off at 
/max = 0.144 arcsec"^ [equation l[9]l]- The full lines are the Han- 
kel transform relations [equation jT}], assuming no low- frequency 
cut-off in the power spectrum. The points and error bars are de- 
rived from 100 realisations of the power spectrum on a 8192^ grid 
with the sampling appropriate for region NPl. The horizontal 
bars represent the bin widths and the crosses the mean separa- 
tion for data included in each of the bins. The error bars represent 
the rms variations for the model structure functions. The dashed 
curves are the structure functions for pure power-law power spec- 
tra with the appropriate low-frequency indices [g = 2.39 and 2.32 
respectively for panels (a) and (b)]. In panel (a), we also show the 
structure function for a power-law power spectrum with g = 11/3 
(dotted line). 
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4.3.2 Modified power-law C{f±) 

In order to fit the observed structure functions, we require 
power spectra of tlie cut-off and broken power-law forms 
given in equations ([SJ and below. For both cases, we 
have verified that the mean of the structure functions from 
realisations of the power spectrum converges to the Hankel 
transform result and that the finite size of the grid for the 
realisations has a negligible effect on the structure functions 
for the separations we sample. Example calculations for both 
types of structure function are shown in Fig. [TT] 



4.4 Observed structure function 

The areas of 3C 31 where we have enough contiguous data 
points to evaluate the RM structure function are discon- 
nected from each other and have very different variances. 
We have therefore analysed the areas NPl, SPl and SP2 
(Fig. |3};) separately at 1.5-arcsec resolution. Use of the 5- 
frequency, 5.5-arcsec RM image (Fig. [2};) allows us to ex- 
amine the larger areas NP2 and SP3, but adds little extra 
information for NPl and SPl, where essentially all of the 
emission is detected at the higher resolution (SP2 and SP3 
overlap). Systematic errors are likely to be important for 
SP2: the signal-to-noise ratio is relatively low and indepen- 
dent RM fitting algorithms give significantly different struc- 
ture functions. We plot the data for this region, but do not 
use them in our determination of the form of the structure 
function. 

In the presence of noise, the structure function has a 
positive bias, and for uncorrelated random noise with rms 
CTflt, its expectation value is 2a^^. We have made a first-order 
correction for this bias by evaluating S{r±) — 5'raw(»'x) — 
2(T|t, where S'raw(''±) is the structure function determined 
directly from the images and am is set to be the rms error 
in the RM for the area in question. In practice, am lies 
between 1.6 and 3radm~^ so the correction is small in the 
North and insignificant in the South. The noise is essentially 
uncorrelated on scales larger than the beam. The measured 
structure functions for the 5 regions are shown in Fig. [I2ja) 
- (e); in Fig.[T2jf), we show a superposition of the structure 
functions for all three southern regions. 

We have investigated the hypothesis that the power 
spectrum (and hence the structure function) has the same 
form for all of the regions, but with varying normalization. 
It is clear from Fig. I12l and from Figs[5Ja) and (b) that there 
is significant power on a range of spatial scales and that a 
Gaussian power spectrum, (for which we expect a flat struc- 
ture function on large scales) is inadequate. Motivated by 
the approximately power-law form for the structure function 
and equation (IBll) . we initially tried an RM power spectrum 
of the form C{f±) oc fj^'' over an infinite frequency range. 
Even after including the effects of the beam [equation (|B2|) . 
Fig. [To], the model structure functions showed insufficient 
curvature to fit our 1.5-arcsec resolution data: all three re- 
gions imaged at 1.5-arcsec resolution require a change in 



3 Note that lEnfilin fc Vogti | |2003|) use wavenumbcr rather than 
spatial frequency and a different Fourier transform convention, 
but their power spectra are the same as ours with the substitution 
= fcx/27r. 



structure-function slope at r± « 10 arcsec; this is also ap- 
parent from the superposition of data for the South of the 
source in Fig. I12f f'). Two functional forms for C{f±) which 
have been discussed in the literature and which fit the data 
reasonably well are a broken power law (hereafter BPL) 

Cif±) = Do.fb""/'(/x//b)"''™ .f± < h 

= DofZ''^^ .f± > ,fb (8) 

and a cut-off power law (hereafter GPL) 

Cif±) = Cofr ,fx</max 

= /x > ,fmax (9) 

The former is motivated by the expectation that the power 
spectrum will have a Kolmogorov form on small scales. 
Fig. [TT] compares the structure functions expected for the 
high- and low-frequency asymptotic forms of these power 
spectra with those corresponding to equations ((S} and ([9]), 
emphasising that there must be a change in slope over the 
observable frequency range. 

For each region, we optimised the normalization over 
a grid of models specified by giow and /b for BPL and q 
and /max for GPL. In order to estimate the errors from 
undersampling, we made multiple realisations of convolved 
RM images corresponding to the BPL model power spec- 
trum on the observed grid and evaluated the rms deviations 
of their structure functions. These represent the errors ex- 
pected purely from statistical fluctuations in the absence 
of noise or systematic errors if the model power spectrum 
is valid everywhere. They are much larger than the errors 
expected from the RM flts, and we therefore used them to 
evaluate chi-squared, which we then summed over all bins 
with separation >FWIIM as a measure of the goodness of 
fit. Note, however, that the errors for adjacent bins are not 
independent. 

In Table[T] we quote the optimised parameters for NP2, 
NPl, SPl and SP3 separately. Uncertainties are approxi- 
mate la (68% confidence) limits for two interesting param- 
eters. Acceptable fits have tightly correlated values of giow 
and /b (or q and /max), so we quote pairs of values defining 
the extrema of contours of constant chi-squared. 

We then determined a best overall flt by minimising chi- 
squared summed over all four regions, giving equal weight to 
each and allowing the normalizations to vary independently. 
We flnd best-fitting parameters of giow = 2.32 and /b = 
0.062 arcsec"^ for the BPL model and q — 2.39 and /max ~ 
0.144 arcsec"^ for the GPL model. The corresponding struc- 
ture functions, including the effects of the convolving beam, 
are plotted in Fig. 1121 where the amplitudes have been scaled 
as in Table [2] to flt the observations (we also include the fit 
for SP2). The rms deviations from multiple realisations are 
plotted as error bars attached to the observed points. 

Fig. [12] shows that the observed structure functions for 
NP2, NPl, SPl and SP3 are in very good agreement with 
the models. Together with the absence of any compelling 
evidence for non-Gaussianity, anisotropy or correlation with 
the small-scale structure of the source (confirmed by visual 
comparison between example realisations, such as those in 
Fig. 1131 and the data) , this indicates that our model power 
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Figure 12. (a) — (e): plots of the RM structure function for the regions indicated in Figs [2l and [3l The horizontal bars represent the bin 
widths and the crosses the mean separation for data included in the bins. The curves are the predictions for the model power spectra 
described in the text, including the effects of the convolving beam. The full and dashed curves denote the broken and cut-off power-law 
models, respectively. The error bars represent the rms variations for structure functions derived using the broken power law power 
spectrum on the observed grid of points for each region. Panel (f ) shows the superposition of structure functions for all three regions in 
the South of the source. 
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Table 1. Structure function parameters for individual regions and for a combined fit to all four regions. The last line of the table refers 
to the combined fit in which the values of qiovj and /b (or g and /max) are the same for all regions, but the individual normalizations are 
varied to minimise the overall chi-squared (these normalizations are listed in Table [2]|- 



Region 


FWHM 






BPL 






arcsec 


Best fit 


Min slope 


Ma 






910W 


/b 


9low 


/b 


9low 


NP2 


5.5 


2.35 


0.178 


2.14 


0.062 


2.54 


NPl 


1.5 


1.84 


0.051 


-0.2 


0.034 


2.45 


SPl 


1.5 


2.76 


0.073 


0.50 


0.024 


3.10 


SP3 


5.5 


2.27 


0.049 


1.70 


0.026 


2.51 


Comb 




2.32 


0.062 


2.15 


0.053 


2.47 



CPL 



slope 


Best fit 


Min slope 


Max slope 


/b 


q 


/max 


9 


/max 


9 


/max 


oo 


2.34 


0.231 


2.10 


0.086 


2.55 


OO 


0.076 


2.35 


0.148 


2.00 


0.122 


2.62 


0.190 


0.155 


2.89 


0.210 


2.46 


0.126 


3.12 


oo 


oo 


2.32 


0.083 


2.03 


0.058 


2.63 


oo 


0.072 


2.39 


0.144 


2.28 


0.134 


2.49 


0.156 



spectra are acceptable representations of the observed RM 
fluctuations. The hypothesis that only the amplitude of the 
power spectrum varies with position is therefore consistent 
with our results, especially given that the approximation 
of constant path length through a region inevitably intro- 
duces some inaccuracies (these are corrected in the 3-d sim- 
ulations of Section [5}. A more rigorous error analysis (best 
done using a Bayesi an maximum likelihood technique, e.g. 
IVogt fc Enfi"iinll2005l ) is outside the scope of this paper. 

The BPL model fits the combined observations very 
slightly better, and we adopt it as our standard, but we can- 
not differentiate between the two functional forms (or sim- 
ilar ones) using our RM data alone. The two model power 
spectra have almost the same form at low spatial frequen- 
cies. The CPL form predicts S{r±) oc at small separa- 
tions (Section IB2|) . quite close to the asymptotic r^^^ ex- 
pected for the Kolmogorov spectrum [equation (|Bip ]. which 
is itself steepened by the effects of convolution. In Fig. 1141 
we show the two power spectra and their associated struc- 
ture functions. The latter are almost identical, particularly 
after convolution (Figs I14t l and f), explaining our inabil- 
ity to distinguish between them. The power in fluctuations 
on scales < 5 arcsec differs significantly for the two models, 
however (Fig. I14b l so in Section 14.51 we attempt to use the 
measured depolarizations at 1.5-arcsec resolution to select 
between them. 



4.5 Constraints from observed depolarization and 
the validity of the short-wavelength 
approximat ion 

Measurement of residual depolarization constrains the am- 
plitude of RM fluctuations on scales smaller than the beam. 
In order to test consistency with our adopted model for the 
power spectrum and to check whether the fitted RM is close 
to the convolution of the true RM distribution with the 
beam [the short- wavelength approximation of equation ((TJ], 
we proceed as follows. 

(i) We make realisations of the RM distributions for 
our adopted power-spectrum parameters on a fine grid 
(0.05 arcsec pixels) and scale the amplitude as appropriate. 

(ii) We then calculate Q and U at each of our observing 
frequencies for a uniformly polarized source. 

(iii) We convolve the resulting images to the observing 



Table 2. Normalizations Dq and Co for the combined fit param- 
eters from Table [T] Observed and calculated depolarization are 
also given for regions NP2, NPl, SPl and SP3. 1: region name. 
2: normalization constant Do for the model power spectrum of 
equation JSjl. 3: normalization constant Cq for the model power 
spectrum of equation (|9]l. 4. Observed Burn law (fc), with errors. 
Only data with p > 4(Tp are included, and for NPl and SPl we 
only include points within 60 arcsec of the nucleus. 5: {k) pre- 
dicted for the BPL power spectrum. 6: As 5, but for the CPL 
power spectrum. Rms sampling errors for the models are also 
given in columns 5 and 6. 



Region 


Do 


Co 




(k) 1 rad^ m 


-4 








Observed 


BPL 


CPL 


NP2 


0.019 


0.61 


46 ±8 


21.8 ±0.7 


23.2 ±0.8 


NPl 


0.064 


1.9 


3±7 


14.7 ±0.1 


9.7 ±0.3 


SPl 


0.26 


7.8 


55 ± 13 


59.7 ±0.9 


40.0 ±1.1 


SP2 


0.25 


7.7 








SP3 


0.61 


19 


372 ± 18 


569 ± 13 


542 ± 7 



resolution and derive the polarized power P = (Cf' + 17^)^^^ 
and E-vector position angle, x- 

(iv) We then fit x against and Inp against A* to make 
images of RM and k as for the observations. 

(v) We average over the images to derive the mean values 
of the Burn parameter (fc) and the degree of polarization 
(p(A)) at each of our observing wavelengths. 

(vi) Finally, we compare the derived RM image with the 
convolution of the full-resolution version with the observing 
beam. 

We have verified that the short-wavelength approxima- 
tion holds to high accuracy for the model power spectra 
in regions NPl, NP2, SPl and SP2 by taking the differ- 
ences between simulated RM images and convolutions of 
the true RM distributions with the appropriate observing 
beams. The rms difference is <0.1radm~^ (negligible com- 
pared with the noise) and the means differ by <0.01 rad m~^ 
in all four regions. Even in SP3, where the short-wavelength 
approximation breaks down, the mean structure function for 
a large number of simulated RM images is in precise agree- 
ment with that derived from the Hankel transform relation. 
The reason is that areas of high depolarization and deviation 
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Figure 13. A comparison between tlie RM data in tlie regions 
used for structure-function analysis and example realisations gen- 
erated using the BPL power spectrum described in the text. The 
mean values for each panel were subtracted before plotting to fa- 
cilitate comparison of smaller-scale structure. The normalizations 
are given in Table[2] (a) and (b): data and model for region NP2. 
(c) and (d); SP3; (e) and (f); NPl; (g) and (h); SPl; (i) and (j): 
SP2. Note that noise in the RM fits is generally uncorrelated on 
small scales, and the observed RM images therefore appear to 
have slightly more small-scale structure than the models. This is 
more obvious in panels (a), (b), (e) and (f), where the amplitude 
of the noise relative to the plotting range is larger. The colour 
scales (±30radm~^ for panels a, b, e and f); ±100radm~^ for 
the remaining plots) are shown by labelled wedges. The horizontal 
lines give a scale of 10 arcsec. 
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Figure 14. (a), (c), (e): the two model RM power spectra dis- 
cussed in the text. Solid line, BPL with indices (Jhigh = 11/3, 5io„ 
= 2.32 and a break frequency of 0.062 arcsec"^ [equation ISllI . 
Dashed line: CPL with q = 2.39 and a high-frequency cut-off at 
/max = 0.144 arcsec"^ [equation JgJ]. (b), (d), (f): structure func- 
tions computed for the power spectra in panels (a), (c) and (d), 
with the same line codes, (a) and (b) no convolution; (c) and (d) 
1.5 arcsec FWHM convolving beam; (e) and (f) 5.5 arcsec FWHM 
convolving beam. 

from rotation occupy a small fraction of the RM image 
and therefore have very little effect on the spatial statis- 
tics (see Appendix [C] for a more detailed analysis of these 
effects) . 

The mean values (fc) predicted for regions NPl and SPl 
by the BPL [equation ((Sjl] and CPL [equation (|9])] models 
for C{f±) are given in Table (2] Both models predict val- 
ues consistent with the observations, so we cannot use the 
measured depolarization to discriminate between them. In 
Figs I15f a) and (b), we compare the observed mean degrees 
of polarization at our six observing wavelengths with the 
predictions of the BPL model derived directly from multiple 
realisations and find very good agreement (the CPL model 
again gives very similar results). 

In regions NP2 and SP3, where we have only low- 
resolution data, the predicted variations of (p) with A are 
not in such good agreement with the observations (Ta- 
ble[2}. For NP2, we predict slightly too much depolarization 
[Fig. llSr c)] . The observed degree of polarization does not fall 
monotonically with increasing wavelength, however, indicat- 
ing either that the approximation of constant intrinsic po- 
larization is invalid or that we have underestimated the sys- 
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Figure 15. Plots of the degree of polarization (on a log scale) 
against A* for the regions: (a) NPl, (b) SPl, (c) NP2 and (d) SP3. 
The filled circles with error bars represent the observed mean de- 
grees of polarization. The crosses show the predicted variations 
of polarization with wavelength for the BPL model with the am- 
plitudes from Table |2] Panel (d) also shows predicted values for 
the BPL models with Dq = 0.38 and /b = 0.062 arcsec"^ (open 
squares) and Dq = 0.61 and /b = 0.037 arcsec"^ (open circles). 
The predicted values for these two cases are almost identical, so 
the symbols appear superposed. In all cases, the intrinsic degree 
of polarization, p(0), has been set to give the best agreement 
(minimum chi-squared) between predicted and observed values. 



tematic errors. In contrast, we predict significantly more de- 
polarization in SP3 than is observed [Fig. I15f dl]i. We show 
in that figure that we can reproduce the run of polarization 
with wavelength by simple modifications to the parameters 
of the BPL power spectrum, for example, by reducing the 
amplitude to Do = 0.38 but keeping the break frequency and 
low- frequency slope at the values in Table [1] or by maintain- 
ing the amplitude and slope at low spatial frequencies but 
decreasing the break frequency to — 0.037 arcsec"^ . These 
discrepancies suggest that the form as well as the ampli- 
tude of the RM power spectrum may vary slightly across the 
source, but we cannot quantify the variation further without 
higher-resolution observations of the RM fluctuations in the 
low-surface-brightness regions NP2 and SP3. 



4.6 The outer scale and autocorrelation lengths 

We have so far implicitly assumed that the outer scale of the 
magnetic-field and RM fiuctuations is large compared with 
the size of the regions over which we evaluate the struc- 
ture function. This is justified by the fact that the struc- 
ture functions continue to rise with increasing separation 



The observed and predicted depolarizations at longer wave- 
lengths in SP3 are less pronounced than p redicted by the Burn 
law, as first pointed out bv lTribbld lll991ar) . 



Table 3. Magnetic autocorrelation lengths for the broken and 
cut-off power-law models [equations JSj and (|9j, respectively]. 
The model used for three-dimensional simulations is the broken 
power law (BPL) with /min = 0.00075 arcsec"^ . 



Afl / kpc 



arcsec ^ 


BPL 


CPL 


0.005 


4.9 


2.8 


0.001 


9.0 


5.7 


0.00075 


10.0 


6.5 



up to rx ~ 100 arcsec. The outer scale cannot be arbitrarily 
large, however, and its value has important implications for 
the estimation of magnetic-field strength and autocorrela- 
tion length. 

Our constraints on the outer scale of the magnetic- 
field fiuctuations from the structure-function analysis are 
weak, because of the poor sampling of large separations. 
For both the broken and cut-off power law models, we infer 
/mill < 0.005 arcsec"^ in regions NP2 and SP3. If the appar- 
ent RM gradient over ~800 arcsec (Fig. |9^) is due to large- 
scale structure in the local field (rather than a gradient in 
Galactic RM) then /min ^ 0.001 arcsec"^ . The autocorrela- 
tion lengths for the magnetic field [As; equation (|A4|| ] and 
RM [Arm; equation (|A5|l ] depend on the outer scale. Equa- 
tions (|A7|) and (|A8|) allow us to determine the autocorre- 
lation lengths for both model power spectra analytically; 
values of As for three representative minimum frequencies 
are given in Table O 



4.7 Origins of the power spectrum 

There is no generally accepted theory of magnetic-field fluc- 
tuations in intra-cluster plasma, although a turbulent cas- 
cade from large to small scales is expected. Dissipation of 
turbulence is predicted to occur on the resistive scale, which 
is far below our resolution limit (e.g. lSchekochihin fc CowlevI 
I2OO6I ) . Energy input might come from sub-cluster mergers or 
galaxy motions within the group and from interactions be- 
tween the radio source and the surrounding plasma. Possible 
input scales L range from L « 25 kpc (the projected sepa- 
ration between NGC 383 and its nearest neighbour and the 
bending radius of the inner jets) to L ~ 500 kpc (the to- 
tal size of the radio source). The magnetic power spectrum 
is expected to have a power-law form with slope q = 11/3 
over some in ertial range if the turbulence is primarily hy- 
drodynamic (|Kolmogorovl Il94l | ) or for the type of MHD 
cascade analysed by Goldreich fc Sridhaij |l993). Our re- 
sults show that the power spectrum could have this form 
on scales <5kpc, but must flatten on larger scales. Such 
a break in slope might occur at a characteristic held rever- 
sal scale Zx ~ 0.01 - O.IL if the field is generated by some 
form of fluctuation dynamo l|Schekochihin fc CowlevI l2006l : 
lEnfilin fc Vogtll2006l ). 
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5 THREE-DIMENSIONAL MODELS 

5.1 General considerations 

We have established that essentially all of the Faraday ro- 
tation associated with 3C 31 is due to foreground material 
(Section [23]). The only known ionized gas component with 
the requisite scale size is the hot plasma associated with 
the su rrounding group of galaxies, which is centred close to 
3C31 ijKomossa fc Bohringeij[l999l ). Although a component 
of Faraday rotation associated with a thin shell of mate- 
rial around the lobes of extragalactic radio sources has been 
postulated by BickncU, Cameron fc Gingold (1990, see also 
iRudnick &i Blundelll |200^, this component cannot domi- 
nate, as there would then be no reason for the path lengths 
to differ between the approaching and receding lobes. There 
would then be no explanation for the systematic asymmetry 
in RM variance and/or depolarization with jet prominence 
seen in 3C31 and many other FRI sources (jMorganti et al.l 
1 19971 ). Cooler (^lO^K) ionized gas has been detected only 
in the inner few arcsec of 3C31 do wen. O'Dea fc Keeilll990l : 
iNoel-Storr et al.ll2003l ) and may be expected to clump on 
scales much smaller than those of the observed large-scale 
RM gradients. We therefore model the Faraday rotation as 
an effect of magnetic field distributed throughout the group 
gas. 

5.2 Simulations with a realistic power spectrum 

Realistic predictions of the RM distribution can be made us- 
in g Monte Carlo simul ations. Our method is essentially that 
of iMurgia et~aLl (|2004l ) : as in our two-dimensional analysis, 
we assume that the magnetic field is an isotropic, Gaus- 
sian random variable. To construct a numerical realisation 
of such a magnetic field in a cubical box, we start from a 
magnetic vecto r potential A(r) with a given power spectrum 
l|Tribblelll991bl ). 

(i) At each point of a 2'" x 2"" x 2"" grid in frequency 
space, we select the real and imaginary parts of the vector 
components of the Fourier transform of the vector poten- 
tial A(f) from a Gaussian random distribution with unit 
variance. 

(ii) We then multiply by the square root of the power 
spectrum of the vector potential. For a magnetic-field power 
spectrum w(/) oc this is oc f~''~^ . 

(iii) Our 2-d analysis was consistent with a magnetic 
power spectrum which has the same form everywhere, but 
varying normalization, so all of the 3-d simulations assume 
the BPL power spectrum of equation (|8]) with giow = 2.32 
and /b = 0.062 arcsec"^. The minimum frequency is set 
to /min ~ 0.00075 arcsec"^, matching the size of the sam- 
pling grids. The magnetic autocorrelation length is As = 
29.6 arcsec (10.0 kpc). 

(iv) The corresponding magnetic field in Fourier space, 
B(f) = 27rif X A(f), is then divergence-free. 

(v) We transform the line-of-sight field component Bz to 
real space and multiply it by the assumed density function 
n(r) to generate a Faraday depth modellf] 

(vi) nBz can then be integrated along the line of sight 



from the emitting surface to the front surface of the cube to 
generate a simulated RM image. 

(vii) Finally, we convolve the RM image with the observ- 
ing beam, assuming that the short-wavelength approxima- 
tion holds, and sample it as for the observations. 

In practice, we used 512^ and 1024'^ grids for the low- and 
high-resolution images, respectively, with cell sizes in real 
space set to 2.6 and 1.3 arcsec. For each of the models de- 
scribed below, we made 25 realisations at low resolution, 
allowing us to explore the variation in predicted profiles, 
together with one realisation at high resolution. We have 
verified that we recover the input power spectrum from 
two-dimensional structure-function analyses for regions over 
which the fluctuation amplitudes are reasonably uniform. 

5.3 Comparison between predicted and observed 
RM distributions 

Our criterion for an acceptable 3-d model is that it fits both 
the form and the amplitude of the RM structure function to 
within errors set by a combination of sampling variance and 
noise in the observations. In order to assess the goodness of 
fit, we need to average over areas which are small enough 
to resolve spatial variations but large enough to give robust 
estimates of the RM fluctuation amplitude. Our approach is 
to make binned proflles of urm at 5.5 arcsec resolution from 
the multiple realisations (sampled in the same way as the ob- 
servations) and then to optimise the overall normalization. 
Our measure of goodness of fit is chi-squared, summed over 
all bins with a significant number of points. As errors, we 
use the rms of the fitting error for the observations (from the 
least-squares analysis of Section I2.3|l and the sampling er- 
ror for the simulations (from the multiple realisations) , com- 
bined in quadrature. We have fit both the 4 and 5-frequency 
RM proflles. The simulations are more tightly constrained 
by the latter, whose fltting errors are significantly smaller. 
Our assumptions are as follows: 

(i) The two sides of the source are intrinsically identical. 

(ii) The emitting and Faraday-rotating media are axisym- 
metric with an ax is orientated at 52.4° t o the line of sight, 
as in the model of lLaing fc Bridl3 (|2002al ) for the inner jets. 

(iii) The surrounding gas density distribution has a radial 
densit y profile as characterised by iKomossa fc Bohringerl 
(| 19991 ). with a central density of 1900 m ^, a core radius 
of Tc = 52 kpc and a form factor /3atm = 0.38. We ignore 
the effects of the gas co mponent associated with the galaxy 
iHardcastle et al"]|2002l . see below). We distinguish two types 
of model: 

Spherically symmetric The source can be represented as 
a plane containing the inner jet axis, i.e. it has negligible 
extent along the line of sight. This is physically unrealistic, 
as we expect thermal matter to be excluded from the volume 
occupied by the relativistic plasma, but should be a good 
approximation where the sour ce is narrow, and is a usefu l 
comparison with earlier work l|Garrington fc Conwaviri99ll : 
iTribblQ,199^ 1. 



^ As noted by iMureia et al" I ||2004| '). this does not precisely pre- 



serve the divergence-free character of the field, but is an adequate 
approximation for our purpose. 
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Cavity Thermal plasma is excluded from the volume oc- 
cupied by re lativistic plasm a, as expected from theoretical 
models (e.g. IScheuerl Il974l ). Cavities apparently devoid of 
X-ray emitting material are indeed observe d around some 
FR I sources (e.g. iMcNamara fc NulsenllioOTl . and references 
therein) We have experimented with various geometries, and 
show the results for ellipsoidal cavities. 

(iv) The normalization of the RM power spectrum then 
depends only on the density n(r) and the rms field B{r). 
We assume that the magnetic field is given by i? oc n^. 
We consider a range from ^ = 0.1 ~ fJ, ~ 1.5. /j, = 0.5 cor- 
responds to equipar tition of magnetic and therm al energy 
and was assumed bv lGarrington fc Conwavl (|l99ll ): /i = 2/3 
is expected for flux-freezing and « 1 is consistent with 
the correlation of X-ray surface brightnes s and RM va riance 
for clusters of galaxies (iDolag et al.il200ll : lDolag||2006l ). The 
best-fitting values of Bo and fj, are correlated for any model, 
since a higher value of fj. gives a relatively lower Faraday 
rotation at large distances from the nucleus, requiring an 
increase in central field strength. 

We show the results for the spherically-symmetric and cavity 
models in Sections 15.41 and 15.51 respectively. 

5.4 Spherically-symmetric distributions of density 
and field 

Figs llGf al and (b) show the comparison between the 
observed profiles of (Jrm and those for the spherically- 
symmetric model. The best fit has a central rms magnetic- 
field strength, Bo = 0.28 nT (2.8 /xG) and fi « 0.5. The fit is 
poor: chi-squared for the 5-frequency profile is 35.7 with 6 
degrees of freedom. Fits of similar quality are obtained for 
any /i in the range 0.1 ^ M ^ 0.8. 

The model profiles show an asymmetry in RM fluctua- 
tion amplitude, but the variation (which occurs on a scale 
^ Tc) is much slower than observed. The predicted RM 
fluctuation amplitude is systematically too high (by a fac- 
tor of two) in the North. It also falls rapidly in the South 
lobe at distances >rc ~ 150arcsec (following the density 
profile), and therefore fails to match the observed value at 
ftiSOO arcsec from the nucleus. 

The failure of this model to fit the remarkably low RM 
fluctuation amplitude in the North of 3C 31 is graphically il- 
lustrated by a comparison between example realisations and 
the observations [Figs llTf a) - (d)]. Corresponding profiles of 
mean RM at two resolutions and rms RM at 1.5 arcsec res- 
olution only are shown in Figs[T8ja) - (f). 

5.5 Cavities 

We have estimated the effect of cavities on the RM vari- 
ation assuming that they are evacuated volumes that do 
not perturb the surrounding gas or field. We first exam- 
ined the effect of " minimal" cavities that match the inner 
jets as modelled bv lLaing fc Bridie (|2002al ). continued out- 
wards. These are cones with half-opening angle 13.2° cen- 
tred on and extending to a distance of 50 arcsec from the 
nucleus, after which they become cylindrical. The effect of 
these "minimal" cavities on the RM distribution in a spher- 
ically symmetric model was too small to be detectable. In 



order for a cavity model to generate the observed rapid RM 
gradient across the nucleus, we found that: (a) the initial 
half-opening angle of the cavities must exceed that of the 
inner radio jet, and (b) the cavities must enclose a volume 
somewhat larger than that of the high-brightness structure 
from which our RM images are derived. 

A simple model which signiflcantly improves the flt 
has ellipsoidal cavities as sketched in Fig. 1191 These have 
semi-major and semi-minor axes of 252 arcsec (85 kpc) 
and 200 arcsec (68 kpc), respectively. The semi-major axis 
projects to 200 arcsec on the plane of the sky. The best flt 
is for jj, ~ 1.0 [Figs I16f c) and (d)]. Example realisations 
and profiles are shown in Figs llTf el and (f) and Figs[l8jg) 
- (i), respectively. The low level of RM fluctuations in the 
North is reproduced quite well, but there is still a signifi- 
cant difference between the observed 5-frequency and sim- 
ulated profiles (chi-squared = 18.6 with 6 degrees of free- 
dom). This is due mostly to one bin, at 300 arcsec from the 
nucleus in the South and coincident with the most nega- 
tive RM values in the source [Figs ^h) and (c)] . The cen- 
tral bins give relatively small contributions to the total chi- 
squared, despite the large differences between observed and 
predicted values, because the sampling errors are large where 
the source is narrow. The 4-frequency profiles of ctrm are for- 
mally consistent with the cavity model for any /i in the range 
0.1</i<1.2 and the chi-squared for /i — 1.0 is 7.0 with 8 
degrees of freedom. 

For this model, the rms central magnetic-field strength 
is Bo = 0.74 nT (7.4 ^G). This is larger than for the 
spherically-symmetric model because there is a lower Fara- 
day depth at large radii for given Bo , partly as a consequence 
of the steeper radial field dependence — 1.0), but mainly 
because there is less gas along the line of sight. Both effects 
require an increase in normalization to match the observed 
Faraday rotation in the South at large radii. The X-ray gas 
parameters were deter mined by fitting to azimutha lly aver- 
aged profiles, however ([Komossa fc Bohringej|l999l ). In the 
presence of a cavity with the shape we assume, the central 
density would be underestimated by a factor ~1.5 and the 
inferred field strength would be overestimated by the same 
factor. 

The predicted ctrm profiles are generally insensitive to 
the precise shape of the cavity provided that it extends to 
the full extent of the sampled region and has a half-opening 
angle comparable with the angle to the line of sight, 9 ~ 52° ; 
for example we have generated similar fits for a cavity which 
initially expands in a cone, becoming cylindrical far from the 
nucleus. The best-fitting value of depends on the precise 
shape of the cavity, since both affect the Faraday depth at 
large radii. 

Although our simple axisymmetric cavity model gives 
a reasonable fit to the observed RM distribution, reality is 
likely to be much more complicated: 

(i) The cavity required to fit the RM data appears much 
larger than the observed extent of the brightest synchrotron 
emission at 1.4 GHz near the nucleus. We note, however, that 
there is low-brightness emission from the North and South 
spurs close to the nucleus and within the projected envelope 
of the cavity. This is visible within 200 arcsec of the nucleus 
on the East side of both jets (Figs[T]and ll9p . 

(ii) It is also possible that earlier epochs of source activity 
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Figure 16. Comparison of observed and simulated profiles of rms Faraday rotation measure ctrm at 5.5 aresee resolution. The blue 
squares represent observed values, with vertical bars corresponding to the rms errors on the RM fit. The red diamonds show the mean 
values from 25 simulated profiles and the vertical bars represent the rms scatter in these profiles. The normalizations of the model profiles 
have been determined by chi-squared minimisation. The curves are the predictions of a single-scale model with D = Xg, as described in 
Section 15.71 The full curve is for an infinite upper integration limit and the dashed line corresponds to truncation of the integral at the 
near surface of the simulation volume, (a) and (b): 3-d model with spherical symmetry (Sec. 15.41 1. The central rms field strength is Bq 
= 0.28 nT and = 0.5. (c) and (d): 3-d model with cavity geometry (Sec. l53t. Bq = 0.74 nT and ^ = 1.0. 



could have cleared a bigger volume of thermal matter. Such 
a cavity would not have been detected in existing Chandra, 
ROSAT or XMM-Newton X-ray images (|Hardcastle et all 
|2002| :Ik omossa fc Bohringej|l999l : ICroston et al.ll2008D . but 
should be clearly visible in a longer XMM-Newton inte- 
gration or, in the radio, at lower observing frequencies (cf. 
iBirzan et a"i]|2008 '). 

(iii) Alternatively, the surprisingly large size required for 
the cavities may just be an artefact of the simplicity of our 
model. Clearly, we cannot constrain the shape of the cavi- 
ties along lines of sight where there is no polarized emission. 
Our inference of a large transverse size for the cavities pro- 
jected on the plane of the sky is purely a consequence of the 
shape of the ctrm profile in those areas we can sample and 



the assumptions of axisymmetry and identical cavities. It is 
possible that the cavities are significantly larger along the 
line of sight than transverse to it or that they differ radi- 
cally from each other in shape. In particular, the structure 
of 3C31 as projected on the plane of the sky is not straight, 
so it is likely that bends also occur in the plane defined by 
the source axis and our line of sight. Such bending might 
well modify the Faraday depth distribution on large scales. 

(iv) If significant magnetoionic material is associated with 
other galaxies in the NGC 383 group located in front of the 
radio emission from 3C31, we might expect to see addi- 
tional, localized, Faraday rotation. An additional contribu- 
tion to the rms RM at a distance of ~300 arcsec from the 
nucleus in the South would greatly improve the (Trm fit, so 
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Figure 17. Comparison of observed and simulated distributions of Faraday rotation measure RM at 5.5 and 1.5 arcsec resolutions, (a) 
and (b): observed, (c) and (d): 3-d model with spherical symmetry (Sec. 15.41 1. (e) and (f): 3-d model with cavity geometry (Sec. 15.51 1. 
The colour scale is the same for all displays but the RM scales for the data are offset by the Galactic RM contribution (Section 12.41 1. 
The bar at the base of panels (a) and (b) shows a 100 arcsec scale. 



22 R.A. Laing et al. 



<RM> 5.5 arcsec 



<RM> 1.5 arcsec 



CTpij 1.5 arcsec 



C 



c 
to 



C 

0) 



-50 



■100 



50 



-50 



50 



-50 



1 1 1 1 1 1 


I I 1 1 11 , 

(a) Data 


, , < 1 < , 






1 1 1 1 1 1 


I I 1 1 11 , 
(d) Spherical- 


1 1 1 1 1 1 


1 1 1 1 1 1 1 




1 1 1 1 1 1 


1 1 1 1 1 1 1 
(g) Cavity 


1 1 < 1 < 1 


1 1 < 1 1 1 1 



I I I I I I I I I I I I I I 

(b) Data 



I I I I I I I I I I I I I 

(e) Sph 



1 1 1 1 1 
- (h) Cav _ 




H 




H 


'in 


H 
H H _ 




H. 


— H- 

H 




H 




H 




H 

1 , , , 


1 



60 



(0 
u 



-a 

(0 



T3 

CO 

S-H 



40 



20 




60 



(c) Ojata 



.III.. I. 



H 1-H 

li.... 



40 



20 - 




60 



(f) Sph 



.III.. I. 



l1i 



40 



b 20 



I ' ' ' I ' 

H : 

(i) Cciv 



-400 -200 200 400 

Distance from nucleus / arcsec 





-200 -200 

Distance from nucleus / arcsec 



Figure 18. Observed and simulated profiles of (RM) and crj^M, with the same binning as in Fig. [9] The observed data at 5.5 arcsec 
resolution are from the 5-frequency fits [Fig. [Sfc)] and the simulated profiles correspond to the example realisations shown in Fig. 1171 
(a) — (c); observed; (d) — (f): 3-d model with spherical symmetry (Sec. [5^ : (g) - (i): 3-d model with cavity geometry (Sec. [STSt . Panels 
(a), (d) and (g) show the mean RM values at 5.5 arcsec resolution; (b), (e) and (h) show mean values at 1.5-arcsec resolution and (c), 
(f) and (i) give the rms RM at the higher resolution. 
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we searched for foreground galaxies close to the RM anomaly 
in the South on the Digital Sky Survey. No candidate objects 
are found within ~2 arcmin, however. 

(v) There could be an entirely "accidental" large-scale in- 
trinsic asymmetry in the foreground Faraday screen that 
produces an excess RM over the outer South jet and tail, 
coincidentally reinforcing the RM asymmetry for the inner 
regions. Some intrinsic asymmetries must be present in the 
gas distribution to provide the observed bending of the radio 
structure, the different shapes of the radio spurs and tails in 
the North and South, and the fla t edge of the East side of 
the South spur (jLaing et al ] |2008l . and this paper, Fig[T} . A 
large-scale asymmetry is also evident in the ROSAT image, 
in which the nucleus of 3C31 is displaced by ~1 arcmin to 
the North-East of the centre of the group gas distribution 
l|Komossa fc Bohringeilll999t ). 

(vi) Finally, our interpretation of the RM profiles close 
to the nucleus might be complicated by Faraday rotation 
ass ociated with the galaxy -scale hot gas component imaged 
bv lHardcastle et al.l l|2002h . which has a core radius of only 
3.6 arcsec. If this had the same field strength and power 
spectrum as we assume for the group gas, then there would 
be a large excess RM close to the nucleus. We have esti- 
mated the resulting ctrm profile and find that there would 
be a detectable excess between —20 and -1-15 arcsec, a factor 
of two increase between —12 and -1-10 arcsec and a peak of 
CRM ~ 400radm~^ within ~1 arcsec of the nucleus (posi- 
tive and negative distances refer to the main (North) and 
counter (South) jets, respectively). We have no information 
about lines of sight within a beamwidth of the nucleus, which 
is unpolarized for our observing frequencies and resolutions 
(perhaps because of such large Faraday rotation gradients 
on small scales). We see no evidence of anomalously high 
RM's close to it, however. In any case, we have no reason 
to suppose that the field strength and power spectrum we 
have adopted are appropriate for the galaxy-scale compo- 
nent, which is smaller than the magnetic autocorrelation 
length we derive for the group gas. We can therefore justify 
our neglect of the galaxy-scale component, noting that it 
may affect the innermost bins of Fig. 1181 

These considerations suggest that it may not be possible - or 
necessary - to explain every feature of the large-scale Fara- 
day RM distribution over 3C 31 with a symmetric model 
of the magnetoionic medium. Nevertheless, we have demon- 
strated that such a model can provide an approximate fit 
to our RM observations and we certainly cannot ascribe the 
entire asymmetry to "accidental" effects and still preserve 
its correlation with jet sidedness in the FRI population as 
a whole. 



5.6 Alternative explanations for the RM 
asymmetry 

Our inference of large cavities surrounding the radio struc- 
ture is testable by deeper radio and X-ray observations. If 
these cavities are not detected, alternative ways to produce 
a Faraday depth distribution consistent with our data must 
be considered. 

One possibility is that the magnetic field is anisotropic. 
A locally two-dimensional structure is expected either if the 
field is "draped" over the surface of a subsonically-expanding 
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Figure 19. The assumed cavity geometry for 3C 31. Top: the 
cavity shape that produced the best fit to the RM data projected 
on a total-intensity image and rotated so the inner jet axis is 
horizontal. Bottom: a cross-section through the cavity in the plane 
defined by the jet axis and the line of sight. The dashed arcs of 
circles represent isodensity contours in the IGM. 

lobe ()Dursi fc Pfrommedl2008l ) or if it is compressed by a 
bow-shock surrounding the source. This could introduce a 
systematic change in Bz/B with position and hence a gra- 
dient of RM fluctuation amplitude. For this mechanism to 
produce a gradient in the observed sense, a special geom- 
etry is needed: the two-dimensional field sheets in front of 
the approaching lobe must be preferentially orthogonal to 
the line of sight while those in front of the receding lobe 
are aligned with it. This might occur if the field is enhanced 
along surfaces that are oblique to the radio jet axis, as in 
a bow shock - such a structure might be imposed at some 
distance from the source by sound/shock waves created by 
earlier outbursts from the AGN (as seen in the Perseus clus- 
ter bv lFabian et al.ll2006l . for example.). 

Alternatively, the Faraday-rotating material might be in 
a disk whose radius (after projection) is similar to, or larger 
than, the sampled area of the tail. A thin disk with a fairly 
uniform distribution of Faraday depth would naturally pro- 
duce a steep gradient of RM fluctuation amplitude across the 
nucleus, coupled with the lack of variation across the outer 
South tail. One possibility is that the gas distribution itself 
is highly flattened. As an example, we have simulated an 
oblate beta-model gas disk with a large core radius contain- 
ing a magnetic field and orientated roughly perpendicular to 
the radio axis. Such a disk could produce the flat proflle of 
CRM in the South, but would also generate higher RM than 
observed in the inner North jet unless it is extremely thin. 
Flattened X-ray emitting gas distributions have occasionally 
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been observed around radio galaxies fe.g. lKraft et al.ir2005h . 
but the group X-ray source surroundins 3C 31 is q uite round 
ijKomossa fc Bohringer|[T999l : ICroston et al.ll2008h . inconsis- 
tent with the appearance of a highly inclined, thin disk. If 
the hot gas is responsible for the Faraday rotation but has 
a spherical distribution, then the only remaining possibil- 
ity is that the magnetic field has a disk-like configuration, a 
supposition that seems entirely ad hoc. 

Finally, the Faraday rotation might come from a dif- 
ferent phase of the ionized IGM than that visualised by 
the X-ray data. Dusty 'superdisks' with diameters > 75 kpc 
and thicknesses Ri25kpc have previously been postulated 
to explain depolarization as ymmetry and other properties 
of powerful radio galaxies (iGopal-Krishna fc Nathl 1 19971 : 
iGopal-Krishna fc WiitalboOOl ). Such a disk in SCSI would 
be too thick to show the rapid change in RM fluctuation 
amplitude across the nucleus that we observe and in any 
case the observed disk of molecular and ionized gas and 
dust is tiny by comparison with the sc ale of the observed 
Faraday rotation - ga2kpc in diameter ( Martel et al.l 1 19991 : 
lOkuda et allbOOSl : iNoel-Storr et allboosi ). This explanation 
for the Faraday-rotation asymmetry also seems extremely 
unlikely. 



5.7 A note on single-scale models 

A rough, but rapidly calculable, approximation to the pro- 
file of rms RM along the source axis can be derived from 
a single-scale mod el of the magnetic-field structure (|Burnl 
ll966l : lFeltenlll996l '). The basic assumption is that the mag- 
netic field is orientated randomly within cells of uniform 
size D and both the strength of the field and the thermal 
electron density are smooth functions of position. The rms 
foreground Faraday rotation (trm, can then be calculated by 
integrating the expression: 



2 _ K^D 

''"RM — 7, 



B {z)n {z)dz 



(10) 



l|Feltenlll996l ) where Zom is the position of the front surface of 
the emitting material along the line of sight z and Zohs is the 
position of the observer (usually Zobs oo). The integral 
does not generally have a solution in closed form, but we 
have verified that our nume rical al gorithm reproduces the 
analytical results derived by iFeltenI (|l99Q) 

in special cases. 

This method allows rapid computation of (trm profiles and 
we have used it extensively to explore different field and 
density combinations. 

For quantitative modelling, the single-scale model has 
serious disadvantages, however. Although it can be used to 
derive a rough value of the field strength, the appropriate cell 
size D to use for a realistic power spectrum is not obvious. 
The apparent RM fluctuation scale, or more precisely the 
RM autocorrelation le ngth Arm [equations (|A5|) and (|A8|) ]. 
is not the correct value. iMurgia et al.l (|2004i l showed that the 
approximation D ~ As, where As is the magnetic autocor- 
relation length [equations (IA4|) and (|A7|l ] holds for power- 
law magnetic power spectra provided that ctrm is measured 
over a sufficiently large area - a restriction which is not sat- 
isfied in our analysis of ctrm proflles at 5.5 arcsec resolution 
close to the centre of the source. We have to limit the size 
of our averaging bins along the jets in order to resolve the 



global variations of path length and RM. Close to the nu- 
cleus, the source structure is dominated by narrow jets, so 
there are only a few autocorrelation areas per bin (Table [3| 
and the single-scale model systematically overestimates the 
rms given by a full 3-d simulation (Fig. I16|l . At distances 
> 100 arcsec, where the source broadens, the single-scale ap- 
proximation with D = Xb gives a much better representa- 
tion of the simulated proflles. Finally, we note that a full 
simulation is required in order to estimate errors for a real- 
istic power spectrum and the sampling grid deflned by the 
observations. 

The single-scale approximation does, however, allow us 
to estimate the effects of the finite depth of the simulation 
volume on the simulated ctrm profile by comparing the pre- 
dictions of equation (|10p with Zobs set to the half-width of 
the simulation box and Zobs = oo. For the cases we have 
modelled in detail, the two profiles are very similar (Fig. ll6|) . 
so almost all of the RM fluctuations will be accounted for 
by the 3-d simulation. 



6 HYDRA A 

Aside from SCSI, Hydra A is the only FRI radio galaxy 
with a well-imaged RM distribution in which there is a 
pronounced asymmetry across the nucleus (|Tavlor fc PerlevI 
[l99S). In SCSI, the orientation is well determined (at least 
near the nucleus) and the group gas is well characterised 
by X-ray observations, but our preferred explanation for the 
observed global RM distribution - that the radio lobes have 
inflated cavities which are devoid of thermal plasma - has 
not yet been adequately tested using X-ray observations. 
For Hydra A, the radio source orientation is less well deter- 
mined, but the cavities associated with th e inner radio lobes 
have been directly ima ged using Chandra (|McNamara et al] 
l2000l : IWise et al.ll2007l '). It is therefore of interest to to ex- 
amine how well a cavity model can account for the global 
RM variations over this source. As for SCSI, we start by 
deriving the RM structure function (Section I6.1|l and use 
the results to simulate the RM variations (Section 16. 2|l . 



6.1 Structure function 

In order to derive an RM power spectrum for our simulations 
and to compare our results with the Baye sian maximum- 
likelihood method of IVogt fc EniSiiiil liooi), we measured 
the RM structure function as in Se ction [H We used RM 
images made using the AlPS R M task (iTavlor fc P crlcv 199- 
and the Pacerman algorithm l|Vogt. Dolag fc Enfilinii200 
as in Section 12.41 We evaluated the structure function for 
the North lobe of Hydra A between 5 and 15 arcsec from the 
nucleus, where both algorithms are in good agreement |f| The 
results derived from t he pacerman RM image are shown in 
Fig. [50] As noted bv lVogt etal] (|2005f ). there is signiflcant 
uncorrelated noise in the RM image, and we have corrected 
the structure function for this using the noise estimates from 
that reference, as in Section 



^ In the South lobe, the RM's are very large, there is significant 
depolarization and the RM distributions calculated using the two 
algorithms differ significantly. We have not, therefore, calculated 
the RM structure function in this part of Hydra A. 
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Figure 20. The RM structure function for the brightest region 
of the North lobe in Hydra A. The curve shows the predicted 
structure function, S{rj_), for a power-law RM power spectrum 
with slope q = 2.77 and a convolving beam of 0.3 arcsec FWHM. 
Error bars show the rms spread for multiple realisations of the 
model power spectrum. 

We fit BPL and CPL models, again as in Section 14.41 
Both fits gave the same solution: a pure power law with 
giow (or q) = 2.77 ± 0.10 with no requirement for a break 
in slope or high-frequency cut-off within our (relatively lim- 
ited) sampling range. The fitted structure function is shown 
in Fig. EOl BPL models with ft > 0.9 arcsec"^ (0.8 kpc"^) 
are a llowed by the data fo r giow ~ 2.7, however. In con- 
trast, IVogt fc Enfiiiiil (|2005l ) found a Kolmogorov spectrum 
for />0.32kpc~^, fiattening at lower frequencies (but they 
did not include the effect of the convolving beam, FWHM 
= 0.32 kpc, in their analysis). 

Thus we find that the index of the power spectrum for 
HydraA (q « 2.8) is steeper than the low-frequency index 
for 3C31 (giow ~ 2.3). Assuming a BPL power spectrum, 
the break frequency is at /{, ~ 0.2kpc~^ for 3C31 compared 
with /f,>0.8kpc~^ for HydraA. 
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Figure 21. Observed and simulat ed RM images for Hy draA. 
(a) Observed RM, using data from lTavlor fc PerlevI ll 19931). The 
RM's for the North and South lobes are from lVogt et al.l ||2005|'| 
and lTavlor fc PerlevI l|l993l ). respectively, (b) A simulation using 
a power-law RM power spectrum with index q = 2.77. The as- 
sumed angle to the line of sight is 45°, B oc n"'^^ (p, = 0.25) 
and the density model includes cavities A and B from lWise et al] 
(2007) as described in the text. The RM scale is from -12500 to 
5000kradm-2. The ima ges have been ro t ated clock wise by 26° 
compa red with those in iTavlor fc PerlevI l|l993) and IVogt et al.l 
ll2005h and the resolution is 0.3 arcsec FWHM. The horizontal 
bars indicate a scale of 20 kpc (18.7 arcsec). 



6.2 Distribution of magnetoionic material 

We attempted to reproduce the large-scale variations over a 
composite RM image of HydraA made using the pacerman 
algorithm for the North lobe, as in Section [6. II and the AlPS 
RM task for the South lob e, wh ere the pacerman solution is 
questionable l|Vogt et al.l 120051 ). There is substantial depo- 
larization in the South, and the AIPS RM image also contains 
a number of artefacts such as large discontinuities on scales 
smaller than a beam. In order to reduce these problems, 
we blanked all points with RM outside the range —15000 to 
+ 10000 rad m~^ and those below a total- intensity threshold. 
The composite RM image is shown in Fig. I21f a). 

We characterised the back ground den s ity di stribution 
using the double beta model of I Wise et al.l (|2003), normal- 



ized to the deprojected profile o f iDavid et ai](j200ll ). The pa- 
rameters of the model are given in Table|4] We ap proximated 
the s hapes of the inner two cavities A and B (|Wise et al.l 
I2OO7I ) as identical ellipsoids with projected semi-major axes 
of 19.8 arcsec exactly along the jet and semi-minor axes of 
11.8 arcsec centred at projected distances of 24.1 arcsec on 
either side of the nuc leus. These are th e mean values for the 
two cavities given bv lWise et all(|2007l ): A and B are on op- 
posite sides of the nucleus, but otherwise very similar in ge- 
ometry and location. Our simulations were as described for 
3C 31 in Section [5.2 1 but taking a power-law magnetic power 
spectrum with q = 2.77, as derived in Section [6. II and an up- 
per frequency limit of /max = 1.67 arcsec"^, set by the size 
of the simulation grid (512^ with 0.3-arcsec pixels). We also 
fixed the lower frequency limit to be /min = 0.012 arcsec"'^. 
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/min is poorly constrained by our structure-function data, 
but the presence of large-scale RM v ariations across the 80 - 
arcsec inner structure of Hydra A (jTavlor fc PerlevI Il993l ) 
requires the power spectrum to have significant amplitude 
at this frequency. 

As the inclination of Hydra A is not well determined, 
we simulated angles to the line of sight in the range 37.5° < 
9 < 60° - larger values failed to generate the observed 
RM asymmetry and smaller ones are improbable given the 
size of the source and its relatively low jet/counter-jet ra- 
tio l|Tavlor et al.lll990l ). We also varied the dependence of 
magnetic field strength on density, taking B oc n'' with 
0.1 < /i < 1.0 For each pair of values 9,ij., we generated 
25 realisations. We then compared the observed and simu- 
lated profiles of ctrm, binning in boxes of length 4.2arcsec 
along the jet axis (PA 26° on the sky) and extending to 
cover the entire width of the source. Finally, we adjusted 
the central rms field strength to optimise the agreement be- 
tween observed and simulated profiles, using a combination 
of fitting and sampling errors to define chi-squared, just as 
for 3C31. 

There is a well-defined minimum in chi-squared, with 
best-fitting values of = 45° and ^ = 0.25. The 68% confi- 
dence region is an ellipse in the 6, jj, plane with semi-major 
axes ae ~ 7° and « 0.2. The agreement between ob- 
served and predicted profiles of ctkm for the best fit is good 
[Fig. I22f a'): chi-squared = 12.7 with 14 degrees of freedom]. 
In particular, the large difference between the fiuctuation 
amplitudes in the two lobes is reproduced well, whereas a 
spherically-symmetric model without cavities gives an un- 
acceptable fit for any inclination. As in 3C31, the discrep- 
ancies close to the nucleus are not significant: errors on the 
predicted values are large because of poor sampling in the 
region of the narrow inner jets. We show an example reali- 
sation in Fig. I21f b). This has been chosen to illustrate that 
the change in sign of RM observed on large scales can be 
generated naturally by our assumed power spectrum with- 
out the need to post ulate an additional ord ered component 
of magnetic field (cf. lTavlor fc PerlevllToQ^ ) . Figs[22l;b) and 
(c) show profiles of observed and predicted mean RM. 

The changes in RM amplitude across Hydra A in this 
picture are due almost entirely to the inner beta model, 
whose core radius is 27.7 kpc (26.5arcsec) and to the cavity 
structure. Any RM fiuctuations due to the larger and more 
tenuous gas component would be essentially constant in am- 
phtude across the central 80 arcsec region which we consider. 
The observed fall-off in RM at distances > 20 arcsec in the 
South lobe is consistent with the drop in density, as pre- 
dicted but not seen in the corresponding (South) lobe of 
3C31. 

The central magnetic field strength required to match 
the RM fiuctuation amplitude for the 45° simulation with 
cavities is 1.9 uT, significantly larger than the 0.7 nT 
deduced from modelling of the North lobe alone by 
IVogt fc EnfilinI (120051 '). We have used the single-scale ap- 
proximation with scale-lengt h Ab to compare our predicted 
profiles of ctrm with those of lVogt fc EnfihnI l|2005h and find 
very good agreement for the North lobe. The discrepancy in 
field strength therefore results from differences in the field 
and density models, as follows: 

(i) The spherically-symmetric model used by 



Table 4. Two-component density model assumed for simulations 
of RM in Hydra A. no is the central density of a component, rc 
is its core radius and /3atm is the form factor. n(r) = no(l + 



no / m 


rc 1 kpc 




6.2 X 10'* 


27.7 


0.686 


2.6 X 102 


235.6 


0.907 



IVogt fc EnfilinI (|2005l) predicts a significantly higher 

density than that of I Wise et"aL I (|2007l ) over the range of 
radii responsible for the bulk of the Faraday rotation. 

(ii) Our inclusion of cavities further reduces the amount 
of gas along the line of sight compared with the expectation 
from a spherically symmetric model. 

(iii) We find that ^ 0.25 gives a better fit to the RM 
fiuctuation pro fi le for a density model with cavities, whereas 
IVogt fc Enfiiiiil (|2005D estimate /i — 0.5 for a spherically 
symmetric model. 

(iv) The power spectrum derived by IVogt fc EnfilinI 
(2005) is significantly different from ours. In order to es- 
timate the effect on the derived central field strength, 
we can comp a re th e magnetic autocorrelation lengths. 
IVogt fc Enfiiiiij (|2005l ) derive As = 1.5 kpc [dividing by 2 
to match our definition from equation (|A4|) ]. compared with 
3.1 kpc for our model power spectrum. 

Our results for Hydra A provide quantitative support 
for our picture of magnetic field distributed through a spher- 
ical background gas distribution with embedded cavities. In 
this source (the only case where the gas density and cav- 
ity structure are both well determined), we can match the 
global variation of RM fluctuation amplitude for plausible 
angles to the line of sight. 



7 SUMMARY AND FURTHER WORK 
7.1 Summary 

We have analysed images of linearly polarized emission 
in the FRI radio galaxy 3C 31 with resolutions of 5.5 
and 1.5 arcsec FWHM at six frequencies between 1365 and 
8440 MHz. Our use of all six frequencies in an analysis of 
the variation of the degree of polarization with wavelength 
has enabled us to measure very small depolarizations, which 
we characterize using the Burn law parameter k. At 1.5- 
arcsec resolution, we find (fc) = 3rad'^m~* (DPicm" = 
0.99) within 60 arcsec of the nucleus in the North jet and 
(fc) = 55rad^m-* iji^f^ = 0.88) in the South. The low 
depolarization and absence of deviation from A^ rotation 
confirm that almost all of the observed Faraday rotation is 
due to foreground material. The RM fiuctuation amplitude 
is significantly larger in the South. 

The amplitude and scale of the RM fluctuations in 
3C 31 appear qualitatively similar to those derived for other 
well-observed FRI sources i n comparable environments , 
e.g. NGC 6251 and 3C 449 (|Perlev. Bridle fc WiUisI Il984l : 
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Figure 22. Profiles of RM along the inner jet axis for observed 
and simulated images of Hydra A (Fig. [2TJ. The simulations as- 
sumed 9 = 45°, Bo = 1.9 nT and /i = 0.25, with the geometry 
described in the text. Data have been averaged in boxes of length 
4.2 arcsec along the jet axis, extending to cover the entire width 
of the source. The jet axis is taken to be vertical in Fig. 1211 or 
PA 26° on the sky. Positive and negative distances refer to the 
North and South lobes, respectively (the former is approaching 
in the simulations), (a) rms RM. The observed values are indi- 
cated by the blue squares, with errors derived from the RM fits. 
The simulation results (red diamonds) represent the mean and 
rms scatter from 25 realisations. The curves show the predicted 
values of ctjjm for a single-scale model with D = \b = 3.1 kpc, 
as calculated for the model power spectrum. Full line: infinite 
integration range; dashed line: integration truncated at the near 
surface of the simulation box. (b) Observed mean RM. (c) Mean 
RM for the simulation shown in Fig. I21f b). This particular real- 
isation (one out of 25) was chosen to illustrate that a large-scale 
change in sign of RM can be generated naturally by an isotropic 
magnetic field with the power spectrum we assume. 



iFeretti et allll999l ). but are much smaller in amplitude than 
those in FRI sources located in cD galaxies with cooling 
cores such as M 87 and Hydra A (Owe n. Eilek fc Kee]|[l99a 
iTavlor fc Perlev, 1993 ) and much larger than for NGC315, 
which is in a sparse group (|Laing et al.ll2006al ). 

We have determined the RM structure function for five 
regions in 3C31 over which the fiuctuation amplitude is 
roughly constant. We model the fiuctuations on the assump- 
tion that the magnetic field responsible for the foreground 
rotation is an isotropic Gaussian random variable (charac- 
terised entirely by its power spectrum). For the range of 
foreground Faraday rotation seen in 3C31, the observed 
RM image is very close to the true image convolved with 
the observing beam. We have used this approximation to 
incorporate the effects of the beam in the calculation of 
theoretical structure functions for model power spectra by 
numerical integration of a Hankel transform relation. We 
evaluate errors due to imperfect sampling by making mul- 
tiple realisations of the model power spectrum on the ob- 
served grid. For the four regions with the most reliable RM 
measurements, our data are consistent with an RM power 
spectrum which has the same form everywhere but whose 
amplitude varies with position. The power spectrum can- 
not be a single power law over the entire range of sampled 
spatial frequencies. A simple functional form that fits our 
data is a broken power law C{f±) oc J^'' with q = 2.32 
for / < 0.062 arcsec"^ and C{f±) oc f]^^^^^, as expected for 
Kolmogorov turbulence, at higher frequencies. A power spec- 
trum with a similar low-frequency slope and an abrupt high- 
frequency cut-off gives almost as good a fit to the observed 
structure functions, however. We can only determine an ap- 
proximate lower limit to the outer scale: in terms of spatial 
frequency, /min ^ 0.005 arcsec"^ (0.015 kpc^^). We conclude 
that these (and other, similar) models of the power spec- 
trum cannot be distinguished using the available RM data 
alone. Our depolarization measurements at 1.5-arcsec reso- 
lution are also consistent with extrapolations of either power 
spectrum to higher spatial frequencies. In the South lobe of 
3C31, the model power spectra (normalized using the RM 
structure function) predict more depolarization than is ob- 
served at low resolution, however. 

In Hydra A, we again fit the RM power spectrum with 
a power law, but measure a steeper slope, q = 2.8, than in 
3C31. We find no evidence for a high-frequency cut-off or 
for a portion of the power spectrum with t he Kolmogorov 
slope at spatial frequencies <0.8kpc~^ (cf. IVogt fc Enfilir] 
l2005h . 

RM analyses for sources in Ab e ll clu sters give qualita- 
tively similar results. iMurgia et all (|2004l ) found a spectral 
slope q ~ 2 extend ing to large spatial scales in Abell 119 
and for Abell 2255, iGovoni et al.l l|2006l ) argued for a spec- 
tral slope which steepens from q « 2 to q « 4 with increasing 
distance from the cluster centre. In contrast. iGuid etti et ah] 
(2008) found that a power-law power spectrum with a Kol- 
mogorov slope and an abrupt long-wavelength cut-off at 
/mill = 0.014 kpc~^ gave a very good fit to their RM and 
depolarization data for A2382, although a shallower slope 
q ~ 2 extending to longer wavelengths was not ruled out. 
All three of these studies assumed high-frequency cut-offs 

/max « 0.1 kpC~^ 

In all cases except Hydra A (where the range of spatial 
scales we have investigated is very small), there is evidence 
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for a change in the power-law slope of the power spectrum, 
but ambiguity about its precise form. Two simple forms of 
the power spectrum which fit all of the datasets are: 

(i) a broken power-law with a high-frequency slope 
Q'high = 11/3 (the Kolmogorov value) and no high-frequency 
cut-off, but flattening to a slope q\o^ « 2 - 3 or even cutting 
off below a frequency /b ~ 0.01 - 1 kpc~^ and 

(ii) a power-law with index g ~ 2 - 3 extending to very 
low frequencies and a high-frequency cut-off at /max ~ 0.1 
- 2kpc-\ 

Other qualitatively similar models would also fit the data. 
Numer ical simulat ions predict a wide range of spectral 
slopes l|Dolaell2006l ). and a detailed comparison between the- 
ory and observation seems premature. 

Our observation of a change in RM fluctuation am- 
plitude across the nucleus of 3C 31 qualitatively supports 
the relativistic-jet interpretation of the brightness asym- 
metry between the jets, wherein the fainter (receding) jet 
is observed through a greater path length of magnetoionic 
medium than the nearer (approaching) jet. To test this 
idea quantitatively, we have modelled the global RM vari- 
ations using three-dimensional simulations, again treating 
the magnetic field as an isotropic Gaussian variable with 
a broken power law power spectrum. We assume that the 
entire structure of 3C 31 has an inclination of 52° as in 
our kinematic models of the inner jets. The profile of the 
RM variance as a function of distance from the nucleus 
of the galaxy is incompatible with a spherically symmetric 
distribution for the magnetoionic medium as derived from 
X-ray observations: this would cause too gradual a varia- 
tion. We expect, however, that the radio source will have 
evacuated cavities in the surrounding hot gas. A model in 
which the edges of the cavities coincide with the bound- 
aries of the radio emission on large scales is more compati- 
ble with the observed RM profile, but current X-ray images 
are not deep enough to reveal any cavities around 3C31. 
The derived central rms magnetic-field strength Bq ~ 0.7 nT 
(7 fiG) and the field varies roughly linearly with density. The 
magnetic field is not dynamically dominant over the vol- 
umes we model. The ratio of thermal to magnetic pressure, 
/3p = Pthermai/(S^/2/io) ~ 10 at the corc radius of the group 
gas, rc = 154arcsec (52kpc). 

We have also modelled the RM distribution in Hydra A, 
where cavities are seen in X-ray images. A three-dimensional 
simulation including the cavities gives a good representation 
of the RM fluctuation profile for an inclination of 45 ± 7 deg 
with B oc fj0.25±o.2 ^ central magnetic field strength of 
1.9nT (19 ^G), so /3p « 30 at = 26arcsec (27.7kpc). 

Close to the nuclei of both sources, our models pre- 
dict higher RM fiuctuation amplitudes than are observed. 
Although the sampling variances are large, this may be an 
indication that our assumption of a power-law dependence 
of magnetic field strength on density is oversimplified and 
that we have therefore overestimated the field strengths over 
these volumes. 

Few other RM profiles have been published for FRI 
sources with jets, and it is not yet clear whether the abrupt 
changes in RM fluctuation amplitude across the nucleus 
observed in 3C 31 and Hydra A are common. In 3C296 
jLaing et al. I l2006bl ), there is a smooth variation of rms RM 
across the nucleus, consistent with the inferred inclination 



of 58° and Faraday rotation produced by a spherically sym- 
metric hot gas component with a large core rad ius, as ob- 
served using XMM-Newton (ICroston e t al."2008Y The RM 
fluctuations observed in NGC 315 iLaing ct al. 2006al) are 
too small to determine a reliable profile, although clearly 
larger on the counter-jet side. 

7.2 Further work 

Our modelling of the RM variations in 3C31 and Hydra A 
suggests a number of observational and theoretical investi- 
gations. 

(i) It is clearly necessary to sample a larger range of spa- 
tial scales to improve constraints on the RM power spec- 
trum and in particular to decide whether it has the form 
expected for Kolmogorov turbulence on small scales. Higher 
resolution would probe small scales directly and remove the 
need to rely on crude estimates from depolarization. Higher 
sensitivity would provide better sampling of large spatial 
scales by allowing study of low-surface-brightness regions. 
Use of analysis techniques such as Bayesian maximum like- 
lihood would give more rigorous error estimates. We have 
shown (Fig. [2]d and c) that it is possible to derive accu- 
rate RM images from a small number of discrete frequencies 
in a suitably chosen observing band: the new generation of 
wideband correlators used by EVLA and eMERLIN will en- 
able simultaneous observations over much larger frequency 
ranges, greatly improving the speed and accuracy of RM 
determination. 

(ii) Large changes in RM fiuctuation amplitude are seen 
across the nuclei in 3C31 and Hydra A, but not in 3C 296. 
The prevalence of this phenomenon and its dependence on 
orientation, large-scale radio structure and environment are 
unknown. It is clear that the RM fiuctuation amplitude 
scales with density, ranging from a few radm~^ in sparse 
environments (e.g. NGC 315) to ~ 10*radm~^ in cooling 
cores such as Hydra A. Observations of a larger sample, em- 
ploying resolutions and frequency ranges matched carefully 
to the environments, will be needed. 

(iii) We have argued that the presence of cavities in the 
external medium must be taken into account when mod- 
elling RM distributions: failure to do so will lead to biased 
estimates of the magnetic field strength and its variation 
with density. Accurate modelling requires knowledge of the 
source geometry (which we can constrain for 3C 31 but not 
Hydra A) and observations of the cavity structure (which are 
available for Hydra A but not yet for 3C31). Deep XMM- 
Newton observations of the 3C 31 group are needed, as are 
comparable data for further sources. Kinematic models of 
the radio jets in Hydra A in other sources could be used to 
constrain their orientations. 

(iv) Deviations in source structure and density from ax- 
isymmetry, Faraday rotation from foreground galaxies, X- 
ray emission from other group members in the field and 
from relativistic electrons deposited in the cavity by the ra- 
dio source (via the inverse Compton process) and variations 
in the Galactic magnetic field across the source structure are 
potential complications which need to be better understood. 

(v) The assumption that the magnetic field is an 
isotropic, Gaussian random variable is questionable. RM im- 
ages of some radio sources show evidence for preferred di- 
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rections and there are theoretical reasons to suppose that 
fields at the edges of cavities and behind radio-source bow 
shocks should be two-dimensional. More detailed simula- 
tions of radio-source evolution in a magnetised intergalactic 
medium are needed. 
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APPENDIX A: NOTATION AND BASIC 
THEORY 

In this Appendix, we summarise the basic theory of RM 
fluctuations in a foreground screen as given in the literatur e 
(Tribble 1991a; Enfilin & Vogt 2003; Mu rgia et al.ll2004l ). 
primarily to establish notation and conventions for compar- 
ison with other work. 



Al Rotation measure 

The rotation Ax of the E- vector position angle of linearly- 
polarized radiation by a foreground magnetized thermal 
plasma is given by; 



Ax = I nB^dz = RMA^ 



(Al) 



where n is the thermal electron density, is the magnetic 
field component along the line of sight, RM is the rotation 
measure and 

K = 



STT^mic^eo 
in SI units. Alternatively, 

Ax/rad = 8.1193 x 10"^ 

X {\/mf / (n/m"^)(B,/nT)d(2/kpc) 



(A2) 



(A3) 



A2 Theory: real space 

In what follows, we make the simplifying assumptions of 
Section [321 Faraday rotation is due entirely to a foreground 
medium in which the magnetic field is a Gaussian isotropic 
random variable and whose density and path length are uni- 
form over the area of interest. 

Our notation is similar to that of lEnfilin fc Vogtlll2003fl . 
r = {x,y, z) is a vector of magnitude r with z along the line 
of sight and rx = {x,y) (magnitude r^) is a vector in the 
plane of the sky. 

The magnetic autocorrelation function wij-) and the 
RM autocorrelation function C(rx) are defined by: 

w{r) = (B(f)-B(f -f-r)) 
C(rx) = (RM(rx-fr'_L)RM(r'_L)) 

and are related by: 

C(rx) = K^n'L i driu{r)r{r^ — r]_)~^^^ 



The autocorrelation lengths for the magnetic field (As) and 
RM (Arm) are of interest for comparison with other work. 
We adopt the definitions: 

A. = Hdr^^ (A4) 
Jo w;(0) 

Akm = [dr,^ (A5) 



These are equivalent to the expression given bv lMurgia et al.l 
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l|2004l ). The values given bv lEnfilin fc Vod (|2003l ) are a fac- 
tor of 2 largerQ We also note that the mean square field 
{B'^) — w(0) and that the average magnetic energy density 
(es) = ui(0)/2/io in SI units. 

A3 Theory: Fourier space 

Here, f = (/^ + fy + fz)^^^ (magnitude /) is a vector in 
frequency space with the coordinate along the line of 
sight. f_L = {fx,fy) (magnitude is a frequency vector in 
the sky plane. We use the Fourier transform convention of 
iBracc wcU (20q3), 

so for one dimension, spatial coordinate 
X and frequency fx the transform pair 17(2;) and g{fx) are 
related by: 



9{fx) 



g(x) exp{~2nifxx)dx 



g{fx) exp{2TTifxx)dfx 



In the isotropic case, we can define the two-dimensional 
RM power spectrum C{f±) such that C{f±)dfxdfy is the 
power for an element dfxdfy in frequency space. The RM 
power spectrum and autocorrelation function then form 
a Hankel transform pair [equation (|6])]. We define the 
three-dimensional magnetic power spectrum w(/) such that 
w{f)dfxdfydfz is the power in a volume dfxdfydfz of fre- 
quency space. This is the Fourier transform of the magnetic 
autocorrelation function and is very simply related to the 
RM power spectrum: 



Mfi-) = 



2g(/x) 



(A6) 



For our model power spectra, it is most straightforward 
to evaluate autocorrelation lengths from integrals over i()(/): 



Arm 



1 /o°°^(/)M 

4/o°°^(/)/^d/ 



(A7) 
(A8) 



The constants in equations (|A7|) and (|A8|) differ by a 
factor of 47r from th ose in the equivalent expressions in 
lEnfiUn fc Voed ()200l their equation 39). 27r of this is be- 
cause we use spatial frequency / rather than wavenumber 
k {k — 27r/) and the remaining factor of 2 is due to the 
definition of autocorrelation length, as explained above. 



trum which has a power-law form C{f±) ~ f^"^ for all fre- 
quencies f±. The structure function then has a simple ana - 
lytical form (jMinter fc Spangle j|l996l : lEnfilin fc Vogtll2003l ). 
In our notation and using equation (|6]): 



2tt 



/i-'Jo(27r/xrx)d/x 



2[C(rx)-C(0)] 

/•OO 

= 4n Jo(2^/xrx)]d/x 
Jo 

Atv"-^ r(2 ^ g/2) 
q~2 r{q/2) 

for 2.5 < g < 4. 

For Gaussian convolution [equation 0]: 



(Bl) 



27r / /l-'Jo(27r/xrx)exp(-2^a71)d/x 



47r / f'-' exp(-2aVl)[l - Jo(27r/x rx)]d/x 
Jo 

2H-<7/2^<j/2^9~2p(2 - q/2) 



q-2 



- 1 



(B2) 



where we have integrated once by parts. We obtain three 
terms: two of these can be integrated using equation 6.631.1 
of Gradshtcyn fc Ryzhik (1994) and combined with equa- 
tion 9.212.2 from the same reference to give a single confiu- 
ent hypergeometric function iFi multiplied by a F function; 
the third is the integral representation of the same F func- 
tion. 



B2 The small-separation limit 

Quite generally, the structure function for a power spectrum 
which is bounded in frequency (unlike the asymptotic case of 
equation IB1|) will be quadratic. S'(rx) oc for suffi ciently 
small separations rx (lTatarskiilll97ll : lTribbl3ll991al ') . For a 
power spectrum C{f±), the structure function is: 



5(rx) 



provided that f±r± <^ 1 over the full range of the integral. 
Similarly for a Gaussian convolving beam. 
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APPENDIX B: ANALYTICAL FORMULAE FOR 
STRUCTURE FUNCTIONS 

Bl Power-law power spectrum 

In order to validate our numerical approach to calculation 
of autocorrelation functions, we first consider a power spec- 



lEnfilin fc Vogtl ll2003h 's definitions of autocorrelation lengths 
are problematic, as they integrate from —00 to +00 in r and r±, 
both of which are intrinsically positive quantities. 



For the specific example of the cut-off power law spectrum 
of equation |5] without convolution, 

^('^^^ " 4^g 

The structure function at small separations then has the 
same functional form as eguation lBll for a power-law index of 
g = 4, leading to possible confusion between power spectra 
with high-frequency cut-offs and steep power laws extending 
to high frequencies. 
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APPENDIX C: ROTATION AND 
DEPOLARIZATION AT LONGER 
WAVELENGTHS 

In this Appendix, we simulate the Faraday rotation and de- 
polarization expected from our model RM power spectrum 
under conditions where the short- wavelength approximation 
does not hold, primarily to aid interpretation of other ob- 
servations. 

In Fig. ICll we show the results of a simulation using the 
method described in Section [4.51 The RM power spectrum 
has the BPL form with the parameters given in Table [1] and 
an amplitude, Do = 0.38, chosen to match the observed de- 
polarization in region SP3. Fig. ICir a) shows the simulated 
RM image convolved to a resolution of 5.5 arcsec. The image 
calculated by evaluating Q and (7 on a fine grid, convolving 
them to 5.5 arcsec resolution, deriving the position angle x 
and fitting for the RM is shown in Fig. lClf bV Although this 
is close to the convolved RM image in most locations, it also 
shows discrepant values associated with large RM gradients 
across the beam. These artefacts are extended perpendicular 
to the direction of maximum RM gradient and are always ac- 
companied by strong depolarization, as shown in Fig. lClf c). 
We see a few artefacts of the predicted type in the observed 
5.5 arcsec RM image, again associated with high RM gradi- 
ents: two examples are shown in Fig. IC2I The reason for the 
artefacts is that there are large deviations from rotation, 
and the RM fitting algorithm does not correct properly for 
the riTT ambiguities in E-vector position angle. This is il- 
lustrated in Fig. IC3I by plots of x against at the points 
labelled a - d in Fig. ICll (these are ordered from least to 
most depolarized). At position a, the RM gradient is small 
and there is negligible deviation from A^ rotation. Positions 
b and c are intermediate cases: Ax is not perfectly propor- 
tional to A^, but the deviations are small enough for the RM 
fitting algorithm to find a reasonable approximation to the 
correct solution. Finally, position d is in a region of high 
RM gradient with large deviations from A^ rotation. The 
fitting algorithm then introduces spurious nvr corrections to 
the position angles and finds an incorrect RM, as is clear 
from Fig. IClf bl . Such deviations from A^ rotation can oc- 
cur even at points with high signal-to-noise ratio, and may 
not, therefo re, be corrected ev en by fitting algorithms like 
PACERMAN (|Dolag et al.l bOOSl ). which take account of the 
spatial coherence of RM. 

For the range of RM, observing frequency and resolution 
appropriate to our observations of 3C31, such artefacts do 
not pose a serious problem, however. Firstly, they should 
cover a small fraction of the area, even in SP3. Secondly, they 
inevitably have low polarized flux density at low frequency 
and are therefore usually excluded from RM fits on grounds 
of low signal-to-noise ratio. 
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Figure CI. Simulated RM and p images generated from a realisa- 
tion of the BPL power spectrum described in the text, computed 
on an 8192 X 8192 grid with 0.05 arcsec pixels. The amplitude is 
Do = 0.38 and the remaining parameters are as given in Table[T] 
The crosses show the locations for the X ^ plots in Fig. lC3l (a) 
Model RM image convolved with a 5.5-arcsec FWHM Gaussian. 

(b) Model RM image of the same area derived by generating Q 
and U images at 5 observing frequencies, convolving them to 5.5- 
arcsec resolution and fitting Ax = RMA^ as for the observations. 

(c) Ratio of the degree of polarization at 1365 MHz to the value 
at infinite frequency, derived from from the convolved Q and U 
images. 
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Figure C2. Part of the observed RM image of the South lobe of 
3C 31 from Fig[2{c). Two examples of artefacts in regions of high 
RM gradient are indicated by white arrows. 




Figure C3. Plots of position angle against for simulated data 
with the BPL RM power spectrum and an amplitude Do = 0.38, 
as required to match the depolarization in region SP3 of 3C31. 
The corresponding RM and depolarization images are shown in 
Fig. ICll where the locations corresponding to individual panels 
are marked. The full lines correspond to the RM's from Fig. lClf a) 
and DP is the ratio of the polarization at 1365 MHz to its intrin- 
sic value. The panels are arranged from least to most depolarized. 
Panel (a) has DP = 0.81 and shows accurate rotation at the 
correct RM. Panels (b) - (d) correspond to positions with higher 
RM gradients, stronger depolarization (DP <C 1) and significant 
deviations from A^ rotation. At position (d), the RM fitting algo- 
rithm has become confused, leading to the RM anomaly visible 
in Fig.EHb). 



